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1.  INTRODUCTION 


The  Interactive  digital  Signal  Analysis  Software  package  (ISAS) 
provides  a  set  of  analysis  programs  which  are  capable  of 
assisting  a  user  in  the  investigation  of  the  dynamic 
characteristics  of  a  system.  The  increasing  use  of  digital 
instrumentation  techniques  and  a  corresponding  increase  in  the 
volume  of  digital  data  have  driven  the  utilization  of  digital 
signal  analysis  software.  The  requirement  to  scan  and  analyze  a 
large  volume  of  data  is  commensurate  with  interactive  analysis 
programs  which  provide  immediate  graphic  representations  of  the 
analysis  results.  A  flexible  methodology  for  the  investigation  of 
a  dynamic  system  is  also  highly  desired.  The  analysis  software 
package  documented  by  this  user's  guide  attempts  to  provide  these 
capabilities  by  supplying  interactive  analysis  tools  that  have  a 
graphic  representation  of  analysis  results  where  appropriate. 
Additionally,  this  user's  guide  exemplifies  a  methodology  for  a 
typical  application  of  the  analysis  software. 


1.1  Background  of  the  Analysis  Software 

The  development  of  the  Interactive  Digital  Signal  Analysis 
Software  Package  resulted  from  a  requirement  to  determine  the 
dynamic  operating  characteristics  of  gas  turbine  engines  during 
propulsion  system  flight  tests.  Digital  data  were  typically 
obtained  at  appropriate  sample  rates  to  capture  the  transient 
event  or  operational  anomaly  of  interest.  At  the  time,  several 
schemes  for  digital  data  analysis  were  available.  However,  the 
associated  analysis  programs  often  imposed  compromising  data 
processing  restrictions  or  required  time  consuming  execution 
through  a  batch  job  processing  structure.  Therefore,  it  was 
deemed  appropriate  to  develop  a  unified  set  of  analysis  routines 
with  interactive  execution  and  graphic  display  of  results. 

The  interactive  digital  signal  analysis  software  package  has  its 
origin  in  the  utilization  of  digital  filter  and  arithmetic 
overlays  for  batch  job  digital  data  processing  available  through 
the  Uniform  Flight  Test  Analysis  System  (UFTAS).  From  lessons 
learned  by  the  application  of  these  overlays  to  propulsion  system 
flight  test  data,  a  foundation  was  established  for  the 
development  of  a  generalized  analysis  methodology  of  dynamic 
data.  This  methodology  consists  of  a  systematic  application  of 
the  interactive  signal  analysis  software  which  is  capable  of 
conditioning,  displaying,  analyzing,  and  filtering  digital  data. 
The  development  of  this  software  consisted  of  the  generation  of 
new  analysis  programs  and  the  enhancement  of  existing  software. 
It  was  necessary  to  develop  methods  of  signal  conditioning  which 
would  correct  data  dropout  and  improper  signal  levels  set  during 
the  digitization  of  analog  signals.  Also,  a  method  to 
interactively  edit  and  plot  time  history  data  for  viewing  on  a 
Tektronix  graphic  terminal  was  developed.  The  standard  UFTAS  file 
structures  and  some  UFTAS  utility  programs  were  employed  for  the 
handling  of  time  history  data.  Many  of  the  frequency  analysis 
programs  were  built  upon  basic  concepts  of  frequency  response 
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analysis  software  algorithms.  Digital  filtering  programs  were 
modified  from  programs  contained  in  a  digital  signal  processing 
software  package  and  enhanced  for  interactive  usage.  The 
resultant  programs  were  consolidated  and  grouped  according  to 
general  function  and  utility.  ISAS  is  maintained  on  a  permanent 
file  fileset  which  is  available  for  execution  on  the  AFFTC  CYBER 
computer  system. 

A  methodology  for  the  utilization  and  interactive  execution  of 
the  analysis  software  is  exemplified  in  this  user's  guide  through 
the  analysis  of  gas  turbine  engine  dynamic  data  obtained  through 
flight  test.  However,  the  analysis  software  is  generalized  and 
can  be  applied  for  the  investigation  of  other  dynamic  systems 
(e.g.  mechanical,  hydraulic,  and  avionics  systems). 


1.2  Example  Application  to  the  Analysis  of  Aircraft  Gas  Turbine 
Engine  Dynamic  Instabilities 

The  complex  design  of  the  modern  gas  turbine  engine  often  harbors 
operational  anomalies.  Operation  of  an  aircraft's  propulsion 
system  at  extreme  regions  of  the  flight  envelope  can  cause 
components  of  the  gas  turbine  engine  to  operate  off  their  design 
point.  Compression  components  can  experience  rotating  stall  or 
surge  which  adversely  affects  other  engine  components.  The 
transient  processes  which  occur  during  a  rapid  advance  of  the 
throttle  may  induce  combustion  instabilities  in  the  engine 
afterburner  and  its  control  system.  Combustion  instabilities 
often  inherent  in  high  performance  thrust  augmentation  systems 
can  result  in  large  amplitude  pressure  oscillations.  These 
pressure  oscillations  can  stall  the  compression  components  of  the 
engine  or  induce  augmentor  blowout.  Phenomena  of  these  types  may 
restrict  the  operation  of  an  aircraft  in  a  region  of  its  flight 
envelope  and  reduce  the  effectiveness  of  its  mission. 

Operational  difficulties  accompanying  engine  component  off-design 
performance  and  throttle  excursions  can  be  minimized  through  an 
understanding  of  the  underlying  physical  principles  and 
characteristics  of  the  encountered  anomaly.  For  example,  an  un¬ 
derstanding  of  the  driving  mechanisms  behind  the  unstable  opera¬ 
tion  of  an  afterburner  can  direct  hardware  modifications  to 
reduce  or  eliminate  the  the  source  of  the  instability.  This 
understanding  of  an  engine  anomaly  is  facilitated  through  the 
performance  of  appropriate  propulsion  system  flight  tests  and  the 
application  of  methodical  post-flight  data  analysis. 

A  methodology  for  the  utilization  of  the  interactive  digital 
signal  analysis  software  to  investigate  unstable  afterburner 
operation  is  exemplified  in  this  user's  guide  through  the 
analysis  of  gas  turbine  engine  dynamic  data  obtained  from  flight 
test.  The  data  analyzed  were  from  a  propulsion  system  which  is 
typical  of  a  modern  high  performance  augmented  turbofan  engine. 
Figure  1  shows  a  cross  sectional  schematic  of  this  type  of 
propulsion  system.  Primary  engine  components  and  stations  are 
indicated.  The  engine  is  typical  of  a  low  bypass  ratio  and  high 


pressure  ratio  design 


Dynamic  engine  instabilities  have  been  observed  to  result  from 
throttle  excursions  which  involve  transient  augmentor  operation 
during  propulsion  system  flight  tests.  These  instabilities  were 
recorded  by  pressure  transducers  located  in  the  engine  core  and 
bypass  flowpaths  and  in  the  the  augmentor  fuel  delivery  system. 
Data  obtained  during  typical  unstable  augmentor  transients  were 
chosen  to  exemplify  the  application  of  the  signal  analysis 
software.  These  augmentor  instabilities  were  characterized  by 
engine  stall  or  stagnation  as  well  as  pressure  oscillations  in 
the  fuel  metering  and  distribution  manifolds.  Analysis  techniques 
typical  of  the  example  application  presented  in  this  user's  guide 
have  revealed  basic  mechanisms  underlying  tnese  unstable 
augmentor  transients.  The  analysis  of  these  gas  turbine  engine 
dynamic  instabilities  was  the  subject  of  an  investigation 
reported  in  Appendix  A  . 


1.3  Structure  of  this  User's  Guide 

This  user's  guide  is  organized  into  several  chapters  which 
explain  the  operational  philosophy  of  the  signal  analysis 
software.  Chapter  2  provides  an  introduction  to  the  analysis 
software  and  explains  the  basic  concepts  underlying  the  structure 
and  control  of  the  individual  programs  and  data  files. 
Functionally  related  software  are  organized  into  program  groups 
which  are  described  in  Chapters  3  through  6.  Each  of  these 
chapters  provides  a  brief  discussion  of  pertinent  technical 
concepts,  a  functional  description  of  the  analysis  routine,  a 
list  of  required  input  and  resultant  output,  a  discussion  and 
example  of  program  execution,  and  a  list  of  limitations  or 
nuances  for  each  analysis  program.  Test  source  data  files  are 
detailed  in  Chapter  7.  These  data  files  primarily  provide  a  check 
for  digital  filter  designs. 


2.  SOFTWARE  STRUCTURE  AND  CONTROL 


The  interactive  digital  signal  analysis  software  is  organized 
into  four  groups  of  programs.  Each  program  group  contains  several 
analysis  programs  which  perform  similar  or  related  tasks  (e.g. 
frequency  response  analysis).  A  detailed  discussion  of  the 
analysis  software  organization,  documentation,  file  standards, 
and  execution  is  given  in  the  following  sections. 


2.1  Fileset  Software  and  Documentation 

The  interactive  signal  analysis  software  package  resides  on  the 
AFFTC  Cyber  170  Model  740  mainframe  computer  system  as  a  direct 
access  NOS  permanent  file  fileset.  This  fileset  is  catalogued 
under  the  username  of  PANDFQ.  The  fileset  ISAS  contains  the 
source,  object,  data,  and  documentation  files  for  the  various 
analysis  programs.  The  program  source  files  are  identified  by  a 
filename  of  the  form  SOURCE/"name",  where  SOURCE  is  the  group 
name  for  the  source  code  of  the  analysis  program  and  "name"  is 
the  element  name  for  the  particular  analysis  program.  Similarly, 
the  executable  object  files  are  identified  by  a  file  name  of  the 
form  ABS/"name"  and  the  data  files  are  identified  by  a  file  name 
of  the  form  DATA/"name".  Associated  with  the  executable  code  of 
each  program  on  the  fileset  is  a  file  of  documentation.  This 
documentation  annotates  the  fileset  directory  listing  and 
provides  a  summary  of  relevant  information  necessary  for  the 
interactive  execution  of  a  program.  This  documentation  is  a  quick 
reference  that  is  available  at  the  terminal  screen. 

The  functional  relationship  of  the  various  analysis  programs, 
associated  data  files,  and  output  files  is  shown  schematical ly  in 
Figure  2.  The  analysis  programs  are  conceptually  organized  into 
four  program  groups  which  are  noted  in  Figure  2.  These  groups  are 
the  data  conditioning  program  group,  the  data  list  and  display 
program  group,  the  frequency  analysis  program  group,  and  the 
digital  filter  program  group. 

The  data  conditioning  program  group  provides  the  following 
capabilities: 

1.  Consolidation  of  parameters  onto  one  data  file  (CONSO) 

2.  Correction  of  instrumentation  and  digitization  errors 
(CONSO  and  SEQCOND ) 

3.  Generation  of  derived  parameters  (SEQCOND) 

4.  Conversion  of  data  files  of  one  format  to  another  (BTOC) 

5.  Temporal  alignment  of  parameter  time  histories  (TIMALIN) 

where  the  name  of  the  associated  analysis  program  is  given  in 
parentheses . 

These  data  conditioning  programs  utilize  the  standard  UFTAS  C- 
file  and  B-file  data  formats  for  the  handling  of  parameter  time 
histories.  This  is  illustrated  in  Figure  2.  These  programs  will 
output  a  resultant  time  history  in  either  a  C-file  or  B-file 
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format.  Utilization  of  these  file  formats  is  discussed  in  Section 
2.2.  The  data  conditioning  program  group  is  discussed  and 
exemplified  in  Chapter  3. 


Another  group  of  related  software  is  the  data  display  and  list 
program  group.  The  following  operations  are  performed  by  this 
program  group: 

1.  List  of  a  parameter  time  history  (LISTF ) 

2.  Plot  of  a  parameter  time  history  (PLOT) 

Figure  2  shows  the  functional  relationship  of  this  software  group 
within  the  software  package.  These  programs  require  that  the 
parameter  time  histories  be  input  in  a  C-file  or  3-file  format. 
Output  from  these  programs  can  be  listed  or  displayed  on  a 
terminal  screen  or  directed  to  a  printer.  This  program  group  is 
discussed  and  exemplified  in  Chapter  4. 

The  frequency  analysis  program  group  provides  the  following 
capabilities: 

1.  Fast  Fourier  transform  of  a  parameter  time  history  (FFT) 

2.  Derivation  of  a  transfer  function  between  parameters 
(TRANS) 

As  Figure  2  indicates,  these  programs  also  require  that  data  be 
input  in  a  B-file  format.  These  programs  produce  results  which 
can  be  graphically  displayed  on  a  Tektronix  graphic  terminal 
screen  or  directed  through  a  local  file  to  be  listed  on  a  line 
printer.  These  programs  are  discussed  and  exemDlified  in  Chapter 
5. 

The  last  program  group  is  the  digital  filter  program  group.  These 
programs  are  shown  in  Figure  2  and  they  provide  the  following 
capabilities: 

1.  Design  of  a  lowpass  filter  file  (LOWPASS) 

2.  Design  of  a  multiple  band  filter  file  (MLTPASS ) 

3.  Implementation  of  a  filter  design  (FILTER) 

The  filter  design  programs  produce  a  filter  file  which  is  based 
on  design  specifications.  If  desired,  a  listing  of  these 
specifications  can  be  displayed  at  the  terminal  screen  or 
directed  to  a  local  file  for  cataloging  or  listing  on  a  printer. 
The  filter  files  are  employed  by  the  filter  implementation  pro¬ 
gram  and  applied  to  data  time  histories  input  in  a  B-file  format. 
This  program  group  is  discussed  and  exemplified  in  Chapter  6. 

The  four  program  groups  provide  analysis  software  capable  of 
processing  time  history  data.  These  programs  and  their  associated 
data,  filter  spec i f ica t i on ,  and  output  files  are  summarized  in 
Figure  2  which  also  presents  a  type  of  flowchart  for  data 
processing.  Paths  of  program  implementation  are  indicated  in  this 
figure.  Many  of  the  processing  paths  are  parallel  and  several  of 
the  analysis  programs  need  not  be  executed  serially.  The  only 


restriction  is  that  data  input  to  an  analysis  program  must  be  in 
either  a  C-file  or  B-file  format.  Beyond  this  requirement, 
implementation  of  the  analysis  software  to  data  time  histories  of 
system  dynamics  is  up  to  the  user's  discretion.  An  example  appli¬ 
cation  of  the  analysis  software  to  observed  gas  turbine  engine 
dynamics  is  presented  in  Chapters  3  thr  „’gh  6.  This  example  will 
illustrate  many  of  the  analysis  paths  through  the  interactive 
signal  analysis  software. 

2.2  Data,  Filter,  and  Output  Files 

Various  types  of  files  are  employed  by  the  analysis  software  to 
communicate  information  to  and  from  the  various  analysis  programs 
within  the  software  package.  The  four  types  of  files  used  are: 

1.  Standard  UFTAS  C-files 

2.  Standard  UFTAS  B-files 

3.  Digital  filter  description  files 

4.  Output  files 

The  relationship  of  these  files  with  respect  to  the  analysis 
programs  is  illustrated  in  Figure  2.  During  processing  of  data, 
these  files  are  attached  to  the  terminal  local  file  space  and  are 
associated  with  the  particular  program  being  executed.  These 
files  can  be  retained  as  permanent  files  and  catalogued  by  the 
NOS  operating  system  of  the  Cyber. 

The  standard  UFTAS  C-files  and  B-files  are  employed  by  the 
analysis  software  to  transfer  time  history  data  from  one  analysis 
program  to  another.  A  B-file  consists  of  a  group  of  parameters 
recorded  over  the  same  time  period  and  written  at  the  same  rate. 
A  C-file  is  a  collection  of  catenated  B-files.  Each  B-file  is  a 
section  of  the  C-file.  C-files  and  B-files  have  unique  header 
records  which  document  a  particular  test  event.  A  detailed 
discussion  of  UFTAS  C-files  and  B-files  is  given  in  Reference  1. 

The  analysis  programs  in  the  ISAS  take  advantage  of  two  key 
features  of  the  UFTAS  B-file  structure.  These  features  are  the 
adjustable  number  of  parameters  available  on  the  file  and  the 
variable  length  of  file.  Figure  2  shows  the  utility  of  using 
these  files  to  communicate  and  store  data. 

Normally,  C-files  and  B-files  are  supplied  to  the  analysis  pro¬ 
grams  by  the  user  or  are  automatically  generated  by  an  analysis 
program  as  the  result  of  data  conditioning.  An  exception  to  this 
convention  are  data  files  supplied  with  the  analysis  programs  in 
the  ISAS.  Currently,  one  such  data  file  exists.  This  data  file 
is  a  B-file  named  DATA/WHITE  and  contains  one  parameter  named 
N0ISE1  which  is  useful  for  testing  a  digital  filter  design.  The 
utility  of  this  program  is  further  explained  in  Section  7.1. 

The  filter  description  files  shown  in  the  right  hand  side  of 
Figure  2  are  used  to  store  design  specifications  of  a  filter. 
These  specifications  are  output  from  the  filter  design  programs 


and  are  employed  by  the  filter  implementation  program  to  filter  a 
parameter  time  history.  These  files  are  similar  to  the  UFTAS 
data  files  and  contain  the  required  number  of  values  for  the 
finite  impulse  response  of  a  digital  filter.  These  files  are 
documented  in  Reference  2. 

Output  files  contain  results  of  an  analysis  program.  These 
sequential  files  are  produced  during  or  at  the  end  of  a  program's 
execution.  They  contain  analysis  results  which  can  be  directed  to 
a  printer  for  listing  or  saved  on  a  permanent  storage  device. 


2.3  Interactive  Execution  and  Documentation 

The  interactive  signal  analysis  software  package  resides  on  the 
AFFTC  Cyber  mainframe  computer  system  as  a  direct  access  fileset. 
This  fileset  is  named  ISAS  and  catalogued  under  the  username  of 
PANDFQ.  The  signal  analysis  software  is  designed  to  operate 
under  the  interactive  facility  of  the  Cyber  NOS  operating  system. 
Some  of  the  analysis  programs  will  provide  graphic 
representations  of  the  analysis  results  as  a  plot  of  the  data 
time  histories  at  a  Tektronix  graphic  terminal.  Otherwise, 
program  execution  and  output  can  be  directed  from  a  nongraphic 
type  of  terminal. 

To  execute  any  analysis  program  at  a  Cyber  remote  terminal,  the 
direct  access  fileset  ISAS  must  first  be  attached  to  the 
terminal  local  file  space  by  entering  the  command: 

ATTACH ,FILESET=ISAS/UN=PANDFQ. 

A  complete  listing  of  the  fileset  elements  (i.e.  source, 
absolute,  and  data  files)  as  well  as  a  one  line  description  of 
the  fileset  element  can  be  obtained  by  entering  the  command: 

LD 

This  command  produces  a  listing  of  the  fileset  elements  at  the 
terminal  screen.  This  listing  is  given  by  Display  1.  A  complete 
description  of  fileset  commands  is  given  in  Reference  3. 

The  executable  code  of  an  analysis  program  can  be  obtained  by 
entering  a  command  of  the  form: 

GF ( ABS/"name" ) 

where  "name”  is  the  particular  name  of  the  desired  analysis 
program.  For  example,  to  access  the  executable  code  for  the  fast 
Fourier  transform  analysis  program  which  is  named  FFT  enter: 

GF (ABS/FFT ) 

FFT  can  now  be  executed  by  a  name  call  load  of  the  program. 
Associated  with  the  executable  code  of  each  analysis  program  is  a 


Display  1  -  Listing  of  the  IDSASP  Fileset  Elements 


file  of  documentation.  This  documentation  annotates  the  fileset 
directory  listing  and  provides  a  summary  of  revelant  information 
necessary  for  the  execution  of  a  particular  program.  This  doc¬ 
umentation  is  a  quick  reference  that  is  available  on  the  terminal 
screen.  The  documentation  for  an  analysis  program  can  be 
displayed  after  the  fileset  has  been  attached  to  the  terminal 
local  file  space  by  entering  a  command  of  the  form: 

GD (ABS/"name” } 

For  example,  the  documentation  associated  with  the  fast  Fourier 
transform  analysis  program  FFT  can  be  displayed  by  entering  the 
command : 

GD (ABS/FFT) 

The  resultant  display  at  the  terminal  screen  is  given  by  Display 
2.  A  complete  listing  of  the  fileset  documentation  is  given  in 
Appendix  B. 
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Display  2  -  Documentation  of  Filesct  K  lenient  ABS/Fl'T 


3.  DATA  CONDITIONING  PROGRAM  GROUP 


The  data  conditioning  program  group  contains  software  which 
consolidates  and  conditions  parameter  time  histories  as  well  as 
manipulates  data  files  in  the  form  of  C-files  and  B-files.  The 
discussion  of  this  software  follows  a  specific  format.  When 
appropriate,  a  brief  technical  discussion  of  the  background  of 
the  analysis  software  is  first  presented.  Program  input 
requirements  are  then  discussed  followed  by  a  discussion  of  the 
interactive  execution  of  the  program.  Next,  typical  program 
output  is  discussed  followed  by  an  example  of  program  execution 
as  applied  to  gas  turbine  pre-stall  data.  Finally,  limitations 
and  nuances  of  the  program  are  given  for  reference.  This  format 
of  discussion  and  example  execution  to  dynamic  engine  parameters 
will  be  repeated  for  each  program  in  Chapters  3  through  6.  It  is 
hoped  that  this  will  provide  confidence  in  applying  the  analysis  < 

software. 


3.1  Data  Circumspection 

Before  application  of  the  analysis  software,  it  is  important  that 
an  all  too  frequently  encountered  truism  "Garbage  In,  Garbage 
Out",  be  considered.  Results  of  an  analysis  can  only  be  as  good 
as  the  data  input  to  the  program.  The  user  of  the  ISAS  should 
be  familiar  with  the  acquisition  and  processing  of  the  data  to  be 
analyzed  from  the  time  it  was  recorded  to  the  time  it  was 
processed  into  an  UFTAS  C-file  or  B-file  format.  Things  to  be 
aware  of  include  aliasing,  time  shifting,  noise,  and  compression. 


Aliasing  is  the  phenomenon  of  high  frequency  energy  being 
transferred  into  lower  frequency  energy  when  analog  signals  are 
sampled  at  improper  rates.  Aliasing  can  have  detrimental  and 
dramatic  effects  on  analysis  results.  To  prevent  this  phenomenon 
from  occurring,  an  anti-aliasing  filter  should  be  placed  before 
the  sampling  device.  An  anti-aliasing  filter  is  an  analog  lowpass 
filter  of  appropriate  cutoff  frequency  such  that  the  higher 
frequency  energy  before  sampling  is  not  aliased  into  the  lower 
frequency  range  of  interest.  A  graphical  approach  to 
understanding  and  preventing  aliasing  is  discussed  in  Reference 
4. 

A  time  shift  of  data  may  occur  due  to  several  sources.  Sensors, 
instrumentation  systems,  anti-aliasing  filters,  and  data  sampling 
may  contribute  to  a  temporal  shift  of  one  signal  relative  to 
another.  Also,  data  processing  by  the  AFFTC  Standard  Analog  to 
Digital  System  (SANDS)  may  result  in  a  temporal  shift  of  digital 
data.  These  effects  can  be  corrected  by  the  analysis  software; 
however,  the  time  shifts  involved  must  be  precisely  known. 

Noise  always  creates  analysis  problems.  Noise  contamination 
should  be  eliminated  or  controlled  at  the  source.  Some  types  of 
noise  may  be  filtered  from  the  contaminated  signal.  This  can  be 
accomplished  by  the  digital  filtering  programs  contained  in  the 
analysis  software  package.  Again,  it  pays  to  closely  follow  data 
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Data  dropout  can  originate  from  many  sources.  The  instrumentation 
system  or  data  handling  operations  may  be  at  fault.  This  phe¬ 
nomenon  can  be  difficult  to  eliminate;  however,  if  the  source  and 
character  of  the  dropout  are  known  appropriate  measures  can  be 
taken  before  processing  the  data.  Data  dropout  has  an  adverse 
affect  on  most  of  the  frequency  analysis  results  and  this  cannot 
be  tolerated.  Data  dropout  is  corrected  by  some  of  the  data 
conditioning  and  handling  programs. 


3.1.1  Example : 

A  discussion  of  the  acquisition  and  processing  of  the  example 
test  data  before  an  analysis  by  the  interactive  signal  analysis 
software  is  given  below.  Also,  the  background  of  the  propulsion 
system  instrumentation  and  post-flight  data  acquisition  system  is 
discussed.  An  understanding  of  these  systems  provides  insight 
into  the  flight  test  data  and  any  anomalies  in  the  data  such  as 
aliasing,  time  shifting,  noise,  or  compression. 

Signals  from  pressure,  temperature,  rotation,  and  position 
sensors  were  recorded  during  propulsion  system  flight  tests.  The 
instrumentation  installed  on  the  augmented  turbofan  engine  is 
shown  in  Figure  3.  The  various  parameters  were  measured  at  the 
locations  indicated  in  Figure  3.  Data  were  recorded  by  a  low 
frequency  response  pulse  code  modulation  (PCM)  system  and  a  high 
frequency  response  frequency  modulation  (FM)  system.  These  data 
streams  were  recorded  on  magnetic  tape  aboard  the  aircraft 
during  the  flight  test.  High  frequency  response  parameters,  such 
as  the  core  and  bypass  duct  flowpath  pressures,  main  fuel  pump 
pressure,  and  the  augmentor  segment  fuel  pressures,  were  recorded 
by  the  FM  system.  Only  the  FM  parameters  are  analyzed  for  the 
examples  of  the  analysis  software.  The  various  FM  parameters  were 
filtered  before  digitization  to  prevent  aliasing  errors. 

Data  from  the  engine  instrumentation  were  processed  by  a  series 
of  digital  computer  programs  which  produced  the  flight  data  file. 
This  process  is  illustrated  in  Figure  4.  During  this  process  the 
FM  parameters  were  digitized  at  a  rate  of  1,000  samples  per 
second.  A  data  base  was  produced  which  was  compatible  with  the 
digital  data  analysis  software.  The  data  base  consisted  of  four 
C-files  which  collectively  contained  ten  parameters.  These  C- 
files  are  represented  by  the  block  of  data  files  at  the  top  of 
Figure  5.  Many  of  the  parameters  analyzed  as  examples  in  the 
following  sections  have  coded  names  that  are  not  indicative  of 
their  origin  or  identity.  A  listing  of  these  parameter  names 
along  with  the  instrumentation  nomenclature  is  given  in  Table  1. 

The  flight  test  results  recorded  on  the  example  C-files 
characterize  the  engine  flowpath  and  fuel  pressure  fluctuations 
experienced  prior  to  engine  stall.  Figure  5  illustrates  the 
methodology  employed  to  analyze  the  pre-stall  test  data.  The 
various  analysis  and  data  transfer  paths  employed  in  the  examples 
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Table  1 


Instrumentation  and  Data  File  Parameter  Nomenclature 


Instrumentation  Parameter  Description 


AJ 

— 

exhaust  nozzle  area  [percent] 

LOD 

— 

augmentor  light-off  detector 

Nl 

— 

fan  rotor  speed  [rpm] 

N2 

— 

core  rotor  speed  [rpm] 

PB 

PP073 

core  burner  inlet  pressure  [psi] 

PF4 

PF4 

core  fuel  pressure  [psi] 

PLA 

— 

power  lever  angle  [degrees] 

PSEG1 

PP293 

augmentor  segment  1  fuel  pressure 

[psi  ] 

PSEG3 

PP295 

augmentor  segment  3  fuel  pressure 

[psi  ] 

PSEG5 

PP297 

augmentor  segment  5  fuel  pressure 

[psi  ] 

PS2.5C 

PP230 

fan  duct  static  pressure  [psi] 

PT2.5C 

— 

fan  duct  total  pressure  [psi] 

PT2 . 5H 

PP227 

fan  core  total  pressure  [psi] 

PS5C 

PP128 

fan  duct  static  pressure  [psi] 

PT6M 

— 

turbine  discharge  total  pressure 

[psi  ] 

VP025 

VP025 

gear  box  vibration  [g] 

VP053  power  shaft  vibration  [g] 


VP058 


are  indicated.  These  paths  will  be  detailed  by  the  discussion  of 
the  example  application  of  the  analysis  software  presented  in 
this  chapter  and  the  following  chapters. 

The  examples  of  the  signal  analysis  software  were  conducted  on  a 
Tektronix  graphic  terminal  to  obtain  a  hard  copy  of  the  graphic 
display  of  the  analysis  results  when  appropriate.  These  displays 
augment  the  example  discussions. 


3.2  Consolidation  of  Parameters  (CONSO) 

CONSO  provides  the  capability  to  consolidate  several  parameter 
time  histories  from  one  to  ten  C-files  into  one  f-file  or  B-file. 
This  capability  is  useful  for  combining  parameter  time  histories 
which  reside  on  many  separate  C-files  due  to  previous  data  pro¬ 
cessing  limitations  or  procedures.  CONSO  will  correct  moderate 
data  dropout  due  to  data  compression. 


3.2.1  Background: 

CONSO  was  developed  as  the  result  of  a  requirement  to  consolidate 
parameter  time  histories  that  were  resident  on  many  separate  C- 
files  which  spanned  equal  time  periods.  This  situation  resulted 
from  the  availability  of  a  limited  number  of  FM  discriminators 
for  the  digititization  of  an  analog  tape.  Also,  CONSO  was  design¬ 
ed  to  correct  data  dropout  due  to  data  compression  which  may  have 
occurred  during  the  generation  of  UFTAS  C-files.  Data  dropout  due 
to  compression  was  prevalent  in  the  engine  test  data  discussed  in 
the  examples  below  and  resulted  when  a  parameter  did  not  change 
in  value  from  one  time  point  to  the  next.  If  possible,  this 
problem  should  be  addressed  during  the  initial  processing  of  data 
into  a  C-file  format.  However,  this  problem  may  persist  and  CONSO 
will  correct  moderate  levels  of  data  dropout. 

The  parameter  merging  philosophy  of  CONSO  is  illustrated  in 
Figure  6.  A  maximum  of  ten  C-files  can  be  input  to  the  program. 
The  parameter  time  histories  to  be  consolidated  must  reside  on 
these  C-files  and  span  equal  time  intervals  at  identical  sample 
rates.  The  parameter  time  histories  can  reside  on  any  B-file 
section  of  the  source  C-files.  A  maximum  of  200  parameters  from 
the  specified  input  B-file  sections  are  then  transferred  from  the 
input  C-files  and  written  over  the  common  time  interval  to  the 
specified  output  B-file  or  C-file.  If  a  C-file  format  is  chosen 
for  the  resultant  time  histories,  the  header  record  from  the 
first  input  C-file  is  written  as  the  header  record  for  the  output 
C-file.  If  data  dropout  is  encountered  during  the  process  of 
transferring  the  parameter  time  histories,  the  previous  value  of 
the  parameter  is  carried  and  used  for  the  dropped  time  point. 
Thus,  execution  of  CONSO  will  produce  a  consolidated  timehistory 
data  file  which  is  free  of  data  dropout  for  subsequent  analysis 
techniques . 


BFEODR  cf 
CFEODR 

C-Flle 


CFLR=C-File  Label 
record 
BFLR=B-File  Label 
record 

BFEODR=B-Fi ' . 

end  of  data  record 
CFEODR=C-File 

end  of  data  record 

ts  =  start  TOD  for 
consolidat ion 
tf-  finish  TOD  for 
consolidation 


B-File 


Figure  6  -  C-File  Consolidation  Methodology  Employed  by  CONSO 


3.2.2  Input: 

CONSO  will  consolidate  parameter  time  histories  which  span 
identical  time  intervals.  The  parameter  time  histories  to  be 
consolidated  must  reside  on  C-files  that  are  attached  to  the 
terminal  local  file  space.  These  files  are  input  to  the  program. 
Specifications  requested  by  the  program's  interactive  dialogue 
are: 

-  number  of  C-files  to  be  consolidated 

-  start  time-of-day  for  the  time  history  interval  of  the 
parameters  to  be  consolidated  [seconds] 

-  finish  time-of-day  for  the  time  history  interval  of 
the  parameters  to  be  consolidated  [seconds] 

-  data  sample  rate  [samples  per  second] 

-  names  of  the  input  C-files 

-  B-file  section  number  which  contains  the  parameters  to  be 
consolidated 

-  selection  of  a  C-file  or  B-file  output  file  format 

-  name  for  the  output  C-file  or  B-file 

Another  consolidation  of  parameter  time  histories  can  be  perform¬ 
ed  upon  consolidation  of  the  specified  C-files.  If  this  option  is 
selected,  the  above  specifications  must  again  be  entered  through 
the  interactive  dialogue.  The  utility  C-file  or  B-file  print 
program  LISTF  may  be  used  to  verify  the  data  consolidation 
performed  by  CONSO. 

3.2.3  Execution: 

Interactive  execution  of  CONSO  is  initiated  by  entering  the 
name  call  load  command: 

CONSO 

The  program  will  respond  by  displaying  at  the  terminal  screen  a 
statement  of  its  capabilities  and  a  prompt  for  an  interactive 
dialogue.  Termination  of  the  interactive  dialogue  returns  the 
terminal  to  the  command  level  of  NOS. 

3.2.4  Output: 

CONSO  produces  one  C-file  or  B-file  which  contains  the 
consolidated  parameter  time  histories.  This  consolidated  file  is 
generated  as  a  local  file  at  the  terminal.  This  file  should  be 
stored  on  a  permanent  storage  device  if  it  is  desired  to  retain 
this  file. 

3.2.5  Example: 

To  facilitate  handling  of  the  engine  pre-stall  test  data,  it  was 
desired  to  consolidate  the  parameters  stored  on  the  four  C-files 


produced  by  SANDS  onto  one  B-file.  This  was  accomplished  by 
utilizing  CONSO  as  indicted  by  Figure  5.  The  parameter  time 
histories  from  the  four  C-files  named  FLT01A,  FLT01B,  FLT01C,  and 
FLT01D  were  consolidated  onto  one  B-file  named  BF001.  This  B-file 
was  subsequently  used  to  list,  plot,  and  edit  the  parameter  time 
histories. 

The  consolidation  of  the  C-file  parameters  is  given  by  a  record 
of  the  interactive  session  with  CONSO  presented  in  Display  3.  The 
parameter  time  histories  from  the  four  C-files  were  consolidated 
over  the  time-of-day  interval  from  75015.000  seconds  to  75018.500 
seconds.  The  four  input  C-files  each  consisted  of  one  B-file 
section.  During  the  consolidation  and  transfer  of  parameters, 
CONSO  corrected  data  dropout.  The  resultant  B-file  named  8F001 
was  generated  as  a  local  file. 


3.2.6  Limitations  and  nuances: 

CONSO  can  consolidate  a  maximum  of  200  parameters.  All  parameters 
on  the  input  C-files  are  consolidated  onto  the  output  C-file  or 
B-file.  Consequently,  if  the  same  parameter  time  history  is 
present  on  more  than  one  input  C-file  it  will  be  repeated  in  the 
resultant  data  file. 
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Display  3  -  Concludod 


3.3  Conversion  of  a  B-file  to  a  C-file  (BTOC) 

BTOC  converts  an  UFTAS  data  file  from  a  B-file  format  to  a  C-file 
format . 


3.3.1  Background: 

This  program  is  a  modified  version  of  the  UFTAS  utility  program 
named  BTOC  which  is  documented  in  Reference  5. 


3.3.2  Input: 

A  local  file  named  BFILE  is  the  primary  input  for  BTOC.  This  file 
must  contain  the  parameter  time  histories  to  be  input  in  a  B-file 
format.  Additiona 1  ly,  a  local  file  named  INPUT  can  be  used  to 
specify  information  to  be  included  on  the  output  C-file  header 
record.  This  information  is  read  by  BTOC  through  the  following 
program  records: 

READ  (  5 / 5 , END=10 )  LIC , ITAI L , ITEST , FLT , DFLT , DREQ , DCOM 
5  FORMAT  (A5,2X,I3,2X,I1,2X,A5,2X,A6,2X,A6,2X,A6) 

READ  (5, ' (10 (A6) ) ' ,END=20)  (T ITLE  ( I ) , I =1 , 10 ) 

The  variable  names  are  documented  in  Reference  5.  Default  values 
for  the  header  record  parameters  are: 

LIC  =  JON13 
ITAIL  *  123 
ITEST  =  1 
FLT  =  FLT13 
DFLT  =  25DEC1 
DREQ  =  25DEC1 
DCOM  =  25DEC1 

TITLE  =  THIS  C-FILE  CAME  FROM  BTOC 

The  program  card  of  BTOC  is  given  below  for  run  time  file 
substitution. 

PROGRAM  BTOC (BFILE ,CFILE , INPUT , OUT PUT , TAPE 1=BF I LE , 
TAPE2=CFILE,TAPE5=INPUT,TAPE6=OUTPUT) 


3.3.3  Execution: 

Interactive  execution  of  BTOC  is  initiated  by  entering  the 
execution  card  statement  as  a  command  at  the  terminal.  This 
name  call  load  command  is: 

BTOC 


Termination  of  BTOC  returns  the  terminal  to  the  command  level  of 


3.3.4  Output: 

Program  messages  are  directed  to  the  terminal  screen. 


3.3.5  Example: 

It  was  desired  to  convert  the  consolidated  engine  stall  time 
histories  from  a  B-file  to  a  C-file  format.  This  format  change  is 
indicated  in  Figure  5.  Data  inspection  through  the  data  list  and 
plot  programs  indicated  that  further  data  editing  as  well  as  the 
derivation  of  a  new  parameter  was  desirable  before  performing 
frequency  analysis  or  digital  filtering  of  the  data. 

Display  4  details  the  interactive  session  of  BTOC  which 
performed  the  data  file  format  conversion.  The  B-file  named  BF031 
was  converted  to  a  C-file  named  CF001.  As  Figure  5  shows,  CF001 
was  employed  for  further  data  processing. 
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Display  4  -  Example  Application  of  the  Data  File  Format  Conversion  Program  BTOi 


3.4  Conditioning  and  Derivation  of  Parameter  Time  Histories 
(SEQCOND) 

SEQCOND  conditions  and  derives  parameter  time  histories  from  a 
standard  UFTAS  time  history  data  file  for  farther  processing  by 
frequency  analysis,  digital  filter,  or  plot  programs.  Several 
data  conditioning  and  manipulation  options  are  available.  These 
include : 

1.  Correction  of  data  dropout  in  a  time  history 

2.  Transfer  of  a  time  history  without  manipulation 

3.  Termination  of  a  time  history  with  a  cosine  taper 

4.  Shift  in  a  time  history  DC  level 

5.  Amplification  or  attenuation  of  a  time  history 

6.  Generation  of  a  derived  parameter  time  history  from  the 
ratio  or  difference  of  two  time  histories 


3.4.1  Background: 

The  data  conditioning  and  manipulation  capabilities  listed  above 
are  implemented  during  the  transfer  of  data  from  an  input  C-file 
to  an  output  B-file.  This  transfer  process  is  interactively 
directed  by  the  user  such  that  only  selected  parameters  of 
interest  over  a  specified  time  period  are  transferred  to  the 
output  B-file.  This  strategy  provides  a  resultant  B-file  which 
contains  only  time  history  data  of  interest  in  a  form  that  is 
compatible  with  the  various  data  analysis  programs  within  the 
program  group.  The  parameter  time  history  conditioning  and 
manipulation  options  of  SEQCOND  are  discussed  below. 

Data  dropout  will  be  corrected  automatically  by  SEQCOND.  The 
program  will  fill  in  dropped  data  points  at  the  appropriate  data 
sample  rate.  The  value  of  the  dropped  data  point  is  set  equal  to 
the  last  valid  data  point.  This  procedure  was  adopted  because 
data  dropout  has  been  found  to  be  primarily  due  to  data 
compression.  Data  dropout  will  be  recovered  if  no  more  than  100 
consecutive  dropped  data  points  are  encountered.  All  parameter 
time  histories  which  are  transferred  to  the  output  B-file  are 
corrected  for  data  dropout. 

SEQCOND  provides  an  option  to  terminate  a  data  sequence  with  a 
cosine  taper  which  is  characteristic  of  the  signal  dynamics 
under  investigation.  This  feature  is  useful  if  "filter  ringing" 
occurs  upon  implementation  of  a  digital  filter  to  a  time  history 
which  has  a  step  or  spike  in  the  signal  level.  Figure  7 
illustrates  the  application  of  a  cosine  taper  to  a  signal  to 
eliminate  a  step  in  the  signal  level.  The  original  signal  is 
retained  until  the  apex  of  the  signal  oscillation  just  prior  to 
signal  dropoff.  Beyond  this  apex  time  point,  the  output  time 
history  consists  of  the  cosine  taper  fit  determined  by  the  period 
and  amplitude  of  the  signal  dynamics. 

It  may  be  necessary  to  shift  the  DC  level  of  a  parameter  time 
history  because  of  instrumentation  or  processing  errors  or  for 
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Figure  7  -  Application  of  Cosine  Taper  Co  Signal  Dropoff 


comparative  analysis  with  another  signal.  Similarly,  a  require¬ 
ment  may  be  to  amplify  or  attenuate  a  time  history  sequence. 
SEQCOND  provides  this  capability  which  is  directed  by  the  program 
user  when  the  data  time  histories  are  transferred  from  a  C-file 
to  a  B-file  format.  These  time  history  conditioning  options  are 
only  applied  to  the  selected  interval  of  the  time  history  for 
output  to  the  B-file. 

SEQCOND  also  provides  a  means  of  deriving  a  parameter  time 
history  from  existing  parameter  time  histories.  Currently,  two 
options  are  available  for  deriving  a  parameter.  These  are  the 
derivation  of  the  ratio  of  two  time  histories  or  the  difference 
between  two  time  histories.  The  derived  parameter  time  histories 
are  appended  to  the  transferred  time  histories  on  the  output  B- 
f  ile. 


3.4.2  Input: 

Parameter  time  histories  to  be  transferred,  manipulated,  or 
derived  must  reside  on  a  single  C-file.  This  C-file  must  be 
attached  to  the  terminal  local  file  space  since  this  data  file  is 
input  into  SEQCOND.  The  program's  interactve  dialogue  is  divided 
into  three  sections.  These  three  sections  require  input  for 
general  processing  of  the  parameters,  input  for  transferring  of 
selected  parameters,  and  input  for  deriving  parameter  time 
histories . 


Specifications  requested  by  the  first  section  of  the  program's 
interactive  dialogue  which  concerns  general  processing  informa¬ 
tion  include: 

-  name  of  the  C-file  which  contains  the  parameter  time 
histories  to  be  processed 

-  section  number  of  the  B-file  which  contains  the  parameter 
time  histories 

-  number  of  parameters  to  be  transferred  or  manipulated 

-  number  of  parameters  to  be  derived 

-  sample  rate  of  the  parameter  time  histories  [samples  per 
second] 

-  start  time-of-day  for  the  output  B-file  [seconds] 

-  finish  time-of-day  for  the  output  B-file  [seconds] 

-  name  for  the  resultant  B-file 

Specifications  requested  by  the  second  section  of  SEQCOND’s 
interactive  dialogue  which  concerns  transferred  parameter 
information  are: 


-  name  of  the  parameters  to  be  transferred 

-  DC  bias  (if  desired)  [parameter  units] 

-  amplification  or  attenuation  factor  (if  desired) 

-  option  of  sequence  termination  with  a  cosine  taper 

If  the  cosine  taper  option  is  selected,  additional  requested 
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specifications  are: 

-  approximate  time-of-day  of  the  signal  step  or  spike  apex 

[seconds]  . 

-  approximate  period  of  the  signal  dynamics  [seconds] 

Specif ications  requested  by  the  third  section  of  the  program's 
interactive  dialogue  which  concerns  derived  parameter  information 
include: 


-  name  for  the  derived  parameter 

-  selection  between  a  ratio  of  parameters  or  a  difference 
between  parameters  to  form  the  derived  parameter  time 
history 

-  reference  parameter  numbers  (these  are  supplied  by  SEQCOND 
in  the  second  section  of  the  interactive  dialogue) 

Upon  conditioning  and  manipulating  specified  parameters  from  the 
input  C-file,  other  time  histories  may  be  conditioned  and 
transferred  to  a  B-file  from  another  C-file.  If  this  is  desired, 
the  above  specifications  must  again  be  entered  into  the  program. 


3.4.3  Execution: 

Interactive  execution  of  SEQCOND  is  initiated  by  entering  the 
name  call  load  command: 

SEQCOND 

The  program  will  respond  with  a  statement  of  its  capabilities  and 
a  prompt  for  an  interactive  dialogue.  Termination  of  SEQCOND 
returns  the  terminal  to  the  command  level  of  NOS. 


3.4.4  Output: 

SEQCOND  produces  a  B-file  which  contains  the  parameter  time 
histories  which  were  transferred,  manipulated,  or  derived  from 
the  parameters  available  on  the  input  C-file. 


3.4.5  Example: 

The  parameter  time  histories  contained  in  C-file  CF001  were 
further  conditioned  and  edited  by  SEQCOND.  Also,  SEQCOND  was 
employed  to  derive  a  parameter  from  the  ratio  of  two  pressure 
parameters  contained  in  CF001.  A  cosine  taper  was  applied  to  the 
transferred  parameter  time  histories  to  eliminate  data  dropoff 
and  subsequently  prevent  filter  ringing.  The  transferred  time 
histories  were  trimmed  in  length  to  1024  points  which  captured 
the  pre-stall  dynamics  of  interest.  As  Figure  5  indicates,  the 
transferred  and  derived  parameters  were  output  on  B-file  BF011 
which  was  used  for  input  to  the  frequency  analysis  and  digital 
filtering  programs. 
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The  execution  and  dialogue  of  SEQCOND  is  exemplified  in  Display  5 
which  is  a  record  of  an  interactive  session.  In  section  1  of 
SEQCOND,  general  specifications  concerning  the  parameters  to  be 
processed  were  given.  These  parameters  were  input  from  section  1 
of  C-file  CF001.  Two  parameters  were  transferred  and  one 
parameter  was  derived  from  these  parameters.  The  sample  rate  of 
the  B-file  parameters  was  specified  as  1000  samples  per  second. 
Also,  start  and  finish  times-of-day  for  the  resultant  time 
histories  were  specified  to  be  75017.001  seconds  and  75018.024 
seconds  respectively.  It  was  desired  to  name  the  output  B-file 
BF011. 

Section  2  of  SEQCOND  required  specifications  of  the  parameters 
that  were  to  be  transferred  or  manipulated  from  the  input  C-file. 
For  each  of  the  transferred  parameters,  PP128  and  PP227,  a  DC 
bias  level  and  amplification  or  attenuation  factor  was  required. 
Parameters  PP128  and  PP227  were  transferred  unchanged  with  a  DC 
bias  of  0  and  an  amplification  factor  of  1  over  the  selected  time 
period.  However,  a  cosine  taper  was  applied  to  these  parameters. 
PP128  and  PP227  were  terminated  with  a  cosine  taper  at  an  apex 
point  of  the  signal  dynamics  which  occurred  at  a  time-of-day  of 
75017.802  seconds.  The  period  of  the  signal  dynamics  at  this 
point  was  approximately  0.019  seconds. 

SEQCOND  serially  assigns  a  parameter  number  to  the  transferred 
time  histories.  This  parameter  number  is  referenced  by  section  3 
of  SEQCOND  which  calculates  derived  parameters. 

Section  3  of  SEQCOND  required  specifications  for  the  parameter 
time  histories  that  were  to  be  derived  from  the  transferred 
parameters.  Parameter  PD001  was  derived  from  the  ratio  of  trans¬ 
ferred  parameter  numbers  1  and  2,  PP128/PP227. 

The  three  resultant  parameter  time  histories  were  written  to  the 
specified  B-file  BF011.  This  B-file  was  output  as  a  member  of  the 
terminal  local  file  space.  The  effect  of  the  cosine  taper  upon 
the  termination  of  parameters  PP128  and  PP227  is  illustrated  by 
plots  of  their  time  histories  which  were  obtained  at  the  terminal 
through  the  plot  program  PLOT. 


3.4.6  Limitations  and  Nuances: 

Data  transfer  from  the  input  C-file  can  include  a  maximum  of  50 
parameters.  From  these  50  parameters  a  maximum  of  25  parameters 
can  be  derived  parameters.  This  limitation  is  system  dependent 
and  may  be  alleviated  in  the  future.  Also,  data  dropout  cannot 
exceed  100  successive  samples  in  a  parameter  time  history. 
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Display  5  -  Example  Application  of  the  Data  Condition  and  Derivation  Program  SEOCOND 


sect  ION  t.  GENERAL  INOU1R*  f  OR  PROCESSING  PARAMETERS 


SECTION  3  INQUIRY  FOR  DEFINED  PAPAHETEFS 


Continued 


PARAMETER  TIME  HISTORY 


Display  5  -  Concluded 


3.5  Temporal  Alignment  of  a  Data  Sequence  (TIMALIN) 

TIMALIN  provides  a  means  of  time  shifting  parameter  time 
histories.  This  capability  is  useful  for  the  temporal  alignment 
of  a  time  history  to  correct  for  time  delays  and  phase  lags 
introduced  by  instrumentation  systems  and  unmatched  filters. 


3.5.1  Background: 

TIMALIN  is  based  upon  a  frequency  response  analysis  (FRA) 
capability  provided  by  the  FRA  program  TSHFT  which  is  documented 
in  Reference  6.  TIMALIN  produces  a  time  history  which  is  lagged 
or  advanced  relative  to  the  input  time  history.  A  time  shift  lead 
is  given  by, 

X  ( 1 0 )  =  x ( t  j  +a ) 

where  x(tp)  is  the  output  time  history  shifted  by  time  a  and 
x  ( t  i )  is  the  input  time  history.  Similarly,  a  time  shift  lag  is 
given  by, 

x  ( 1 0 )  =  x(tj-a) 

The  temporal  shift  of  the  input  time  history  is  accomplished  by 
first  applying  a  fast  Fourier  transform  to  the  time  history.  The 
resulting  Fourier  coefficients  are  then  phase  shifted  according 
to  the  desired  time  shift.  An  inverse  fast  Fourier  transform  is 
then  applied  to  the  shifted  Fourier  coefficients  to  restore  a 
time  history  of  shifted  data. 

The  input  time  history  must  be  bandlimited.  Also,  TIMALIN 
requires  that  the  input  time  history  consist  of  a  number  of 
points  which  is  a  power  of  two.  TIMALIN  will  automatically  zero 
fill  an  input  time  history  to  the  next  power  of  two  if  this  re¬ 
quirement  is  not  met. 


3.5.2  Input: 

TIMALIN  requires  that  the  time  history  to  be  temporally  shifted 
reside  on  a  B-file  at  the  terminal  local  file  space. 
Specifications  requested  by  TIMALIN  through  the  interactive 
dialogue  are: 

-  name  of  B-file  which  contains  the  parameter  time  history 
to  be  time  shifted 

-  name  of  the  parameter  to  be  time  shifted 

-  start  time-of-day  of  the  time  history  [seconds] 

-  finish  time-of-day  of  the  time  history  [seconds] 

-  sample  rate  of  the  time  history  [samples  per  second] 

-  desired  time  shift  (positive  for  lead,  negative  for  lag) 
[seconds] 

-  name  of  the  output  B-file  to  contain  the  shifted  parameter 

-  option  to  write  the  shifted  parameter  as  a  unique 


parameter  or  to  overwrite  the  input  parameter  with  the 
shifted  parameter 


If  the  shifted  parameter  is  to  be  output  as  an  unique  parameter, 
additional  requested  information  is: 

-  name  of  the  shifted  parameter  time  history 

Upon  time  shifting  one  parameter,  if  it  is  desired  to  shift 
another  parameter,  the  above  specifications  must  be  entered 
again. 

3.5.3  Execution: 

Interactive  execution  of  TIMALIN  is  initiated  by  entering  the 
name  call  load  command: 

TIMALIN 

The  program  will  respond  with  a  statement  of  its  capabilities  and 
a  prompt  for  an  interactive  dialogue.  Termination  of  TIMALIN 
returns  the  terminal  to  the  command  level  of  NOS. 


3.5.4  Output: 

TIMALIN  produces  a  local  B-file  which  contains  the  shifted  time 
history  as  a  unique  parameter. 

3.5.5  Example: 

After  the  digital  filtering  of  parameter  time  histories,  it  may 
be  desirable  to  correct  a  time  delay  in  the  time  histories  due 
to  the  linear  phase  response  of  the  filter.  As  Figure  5 
indicates,  the  parameter  time  histories  filtered  by  FILTER 
through  the  MLTPASS  filter  design  were  temporally  realigned  by 
TIMALIN.  This  restored  the  filtered  parameter  time  histories  into 
alignment  with  the  unfiltered  parameters.  This  facilitated 
analysis  of  the  time  histories  in  the  time  domain. 

The  execution  and  dialogue  of  TIMALIN  is  exemplified  in  Display  6 
through  its  application  to  the  filtered  parameter  FP128M 
contained  on  the  B-file  named  BF024.  FP128M  was  shifted  by  a  time 
lag  of  0.032  seconds  over  the  interval  time-of-day  from  75017.001 
seconds  to  75018.024  seconds.  As  Sections  6.2  and  6.3  discuss, 
the  application  of  the  MLTPASS  filter  design  to  the  parameter 
PP128  resulted  in  a  time  history  FP128M  that  was  shifted  forward 
in  time  by  0.032  seconds  relative  to  the  input  time  history.  The 
corrected  time  history  was  output  as  a  unique  parameter  named 
SP128M  and  written  to  the  local  B-file  BF031.  The  effect  of 
TIMALIN  is  shown  in  the  plots  of  parameters  FP128M  and  SP128M 
which  are  presented  in  Display  6. 
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TIME  OF  DAY  [SECONDS! 


PARAMETER  TIME  HISTORY 


4  DATA  LIST  AMD  DISPLAY  PROGRAM  GROUP 

The  data  list  and  display  program  group  contains  software  which 
lists  and  plots  parameter  time  histories.  As  Figure  2  indicates, 
these  programs  obtain  parameter  time  histories  from  UFTAS  C-files 
and  B-files.  These  programs  also  produce  output  in  the  form  of  a 
graphical  display  at  a  Tektronix  graphic  terminal  or  a  listing 
which  can  be  displayed  at  the  terminal  screen  and  written  to  a 
local  file  for  printing  at  a  printer. 

4.1  List  of  a  Parameter  Time  History  (LISTF) 

LISTF  is  a  C-file  or  B-file  time  history  print  program.  LISTF 
will  write  a  listing  of  the  input  file  to  the  terminal  screen  or 
to  a  local  file  for  printing  at  a  remote  site. 

4.1.1  Background: 

LISTF  is  a  modified  version  of  an  CJFTAS  utility  program  named 
LISTBC.  LISTBC  is  documented  in  Reference  4.  LISTF  provides 
options  which  enable  the  specification  of  the  data  print  rate, 
the  parameters  to  be  printed,  and  the  format  of  the  printed  data. 
These  options  are  specified  through  an  execution  command.  This 
command  enables  interactive  execution  of  the  program.  Batch 
execution  of  LISTF  is  possible  and  is  also  documented  in 
Reference  4. 


4.1.2  Input: 

Input  to  LISTF  is  provided  through  an  execution  command  which  is 
entered  at  the  terminal  and  begins  interactive  execution  of  the 
program.  The  execution  command  is  of  the  form: 

LISTF  (FILE,PR=nl,PC=n2,SD=n3) 

where, 

FILE  =  name  of  the  local  file  which  contains  the  parameters  to  be 
listed.  Default  is  the  local  file  named  CFILE. 

PR  =  print  rate  which  specifies  every  nlth  point  to  be  printed. 
Default  is  nl=l  which  prints  every  data  point. 

PC  =  flag  to  indicate  that  only  specified  parameters  are  to  be 
printed.  If  n2  =  0  (default)  all  the  parameters  on  the 
specified  B-file  or  C-file  are  printed.  If  n2=l  the  names 
of  parameters  to  be  printed  are  read  from  a  local  file 
named  INPUT.  The  names  of  the  parameters  are  specified  on 
this  file  in  free  format  in  columns  1  through  80  with  a 
blank  or  comma  as  terminators.  An  example  for  file  INPUT 
is : 
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*EOR 
PI  P2 
P3,P4 
P5 

*EOR 

This  results  in  the  parameters  PI,  P2,  P3,  P4,  and  P5 
being  printed,  provided  that  they  exist  on  the  input  local 
file. 

SD  *  flag  to  indicate  the  number  of  significant  digits  for  the 
printed  data  values.  If  n3=0  (default),  5  significant 
digits  are  printed. 

The  arguments  of  the  LISTF  execution  command  are  optional. 


4.1.3  Execution: 

Interactive  execution  of  LISTF  is  initiated  by  entering  the 
execution  command  at  the  terminal.  If  execution  defaults  are 
invoked,  then  list  is  executed  by  entering  the  name  call  load 
command : 

LISTF 

A  listing  of  the  parameter  time  histories  will  be  directed  to  the 
terminal  screen.  The  program  card  of  LISTF  is  given  below  for  run 
time  file  substitution. 

PROGRAM  LISTF  (CFILE ,TAPE1=CFILE , INPUT , OUTPUT , 

TAPE5=INPUT ,TAPE6=0UTPUT , TAPE 99 ) 

Termination  of  LISTF  returns  the  terminal  to  the  command  level  of 
NOS. 


4.1.4  Output: 

If  run  time 
1 isting  will 


file  substitution  is  not  employed, 
be  directed  to  the  terminal  screen. 


formatted 


4.1.5  Example: 

As  Figure  4  illustrates,  UFTAS  produced  a  data  base  compatible 
with  the  digital  signal  analysis  software.  This  data  base 
consists  of  C-files.  A  listing  of  the  parameter  time  histories 
obtained  from  the  engine  pre-stall  test  data  using  LISTF  was 
desired  for  a  cursory  inspection  of  the  data.  Data  processing 
anomalies  such  as  data  dropout  were  detected  from  these  listings. 
Also,  at  many  intermediate  steps  in  the  example  data  analysis  it 
was  advantageous  to  obtain  a  listing  of  a  parameter  time  history 
for  inspection  at  the  terminal  screen. 


Utilization  of  LISTF  for  the  pre-stall  data  is  indicated  in 
Figure  5.  List  was  employed  to  list  the  data  base  C-files  as 
shown  by  Display  7.  A  listing  of  the  parameters  PP230  and  PP127 
contained  on  C-file  FLT01B  was  obtained  at  the  terminal  screen  at 
the  specified  print  rate  of  100  (i.e.  print  every  100th  point). 
LISTF  was  also  employed  to  confirm  the  consolidation  of 
parameters  performed  by  CONSO  as  discussed  in  Section  3.2. 
Display  8  is  a  record  of  the  interactive  usage  of  LISTF  to 
produce  a  listing  of  the  consolidated  B-file  BF001  at  the 
terminal  screen.  The  ten  consolidated  parameters  are  listed  at  a 
print  rate  of  1,000. 
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Display  7  -  Example  Application  of  the  Parameter  Listing  Program  LISTF  to  a  C-File 
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Display  7  -  Concluded 


Display  8  -  Example  Application  of  the  Parameter  Listing  Program  LIS1F  to  a  B-File 


Display  8  -  Concluded 


4.2  Plot  of  a  Parameter  Time  History  (PLOT) 

PLOT  will  display  a  plot  of  a  B-file  parameter  time  history  on 
the  screen  of  a  Tektronix  graphic  terminal. 

4.2.1  Background: 

PLOT  provides  a  plot  of  a  B-file  time  history  through  interactive 
commands  entered  at  a  Tektronix  graphic  terminal.  The  user  can 
dynamically  modify  the  resultant  plot.  This  capability  provides  a 
means  to  examine  in  detail  many  aspects  of  signal  dynamics  and 
to  estimate  character istics  of  the  dynamics  such  as  the  amplitude 
and  period  of  oscillations. 

4.2.2  Input: 

The  parameter  time  histories  to  be  plotted  must  reside  on  B-files 
which  are  attached  to  the  terminal  local  file  space.  These  files 
are  input  to  PLOT.  Specifications  required  by  the  program's 
interactive  dialogue  include: 

-  model  number  of  the  Tektronix  graphic  terminal 

-  name  of  the  B-file  which  contains  the  time  history  to  be 
plotted 

-  name  of  the  parameter  to  be  plotted 

-  units  of  the  parameter 

-  data  sample  rate  [samples  per  second] 

-  plot  start  time-of-day  [seconds] 

-  plot  finish  time-of-day  [seconds] 

-  option  to  print  a  line  of  remark  at  the  top  of  the  plot 

-  option  to  invoke  automatic  scaling  of  the  plot  ordinate 

If  manual  scaling  of  the  plot  was  selected,  additional  required 
specifications  are: 

-  maximum  value  of  the  plot  ordinate 

-  minimum  value  of  the  plot  ordinate 

The  maximum  and  minimum  value  of  the  plot  abscissa  are  determined 
by  the  plot  start  and  finish  tiraes-of-day. 


4.2.3  Execution: 

Interactive  execution 
name  call  load  command: 


PLOT  is  initiated  by  entering 


The  program  will  respond  with  a  statement  of  its  capabilities  and 
a  prompt  for  an  interactive  dialogue.  Termination  of  PLOT  will 
return  the  terminal  to  the  command  level  of  NOS. 
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4.2.4  Output: 

PLOT  produces  a  plot  of  a  parameter  time  history  on  the  screen  of 
a  Tektronix  graphic  terminal.  A  copy  of  the  plot  can  be  obtained 
through  the  usage  of  a  Tektronix  hard  copy  unit. 


4.2.5  Example: 

Plots  of  the  FM  parameter  time  histories  revealed  details  of 
augmentor  transient  phenomena  which  preceded  engine  stall.  An 
understanding  of  the  onset  of  pressure  fluctuations  was 
facilitated  by  the  high  digitization  rate  of  1000  samples  per 
second . 

The  plot  of  two  B-file  parameters  is  exemplified  by  an 
interactive  session  of  PLOT  as  given  by  Display  9.  The 
parameters  PP128  and  PP227  were  input  from  the  B-file  named  BF001 
and  plotted  on  the  screen  of  a  Tektronix  graphic  terminal.  This 
is  indicated  in  Figure  5.  The  parameter  time  histories  were 
plotted  over  a  time-of-day  interval  from  75015.000  seconds  to 
75018.500  seconds.  Automatic  scaling  was  chosen  for  parameter 
PP128  while  the  plot  ordinate  of  parameter  PP227  was  manually 
scaled  from  0  psi  to  25  psi. 

The  plots  of  Display  9  show  the  onset  of  pressure  fluctuations  in 
the  fan  duct  which  immediately  preceded  engine  stall  at 
approximately  75018  seconds.  During  this  period  of  time, 
augmentor  fuel  system  parameters  indicated  stabilized  fuel  flow. 
Also  duct  flow  pressure  oscillations  were  observed  at  the  high 
pressure  compressor  inlet  but  not  at  the  exit  of  the  high 
pressure  compressor.  This  information  suggested  that  the  large 
amplitude  pressure  oscillations  may  be  associated  with  a 
combustion  instability  of  the  augmentor.  This  situation  is 
summarized  by  a  compilation  of  time  histories  presented  in  Figure 
8.  The  elapsed  time  interval  in  Figure  8  corresponds  to  the 
times-of-day  in  Display  9  from  75017.200  seconds  to  75018.200 
seconds . 

The  augmentor  transient  time  histories  have  suggested  several 
phenomena  which  involve  oscillations  of  the  augmentor  fuel 
delivery  system  and  induced  oscillations  of  the  fan  and  core 
flowpath.  These  phenomena  appear  very  similar,  occur  within  short 
time  intervals,  and  result  in  different  operational  anomalies.  It 
was  not  certain  by  viewing  the  data  time  histories  if  the 
pressure  oscillations  in  the  engine  flowpaths  resulted  from  the 
same  phenomena,  coupled  phenomena,  or  separate  mechanisms.  A 
quantitative  method  for  analyzing  the  time  history  data  was 
needed.  One  approach  was  to  analyze  the  time  history  data  in  the 
frequency  domain.  This  approach  is  exemplified  in  Chapter  5. 
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Figure  8  -  Time  History  of  FM  Data  Preceding  Engine  Stall 


THIS  INTERACT  I VE  P ROC RAH  WILL  PLOT  A  USER  SPECIFIED  PARAMETER 
SEOUENCE  VERSUS  TIME  THIS  SEOUENCE  MUST  BE  JN  A  B-FILE  FORMAT. 
PLOT  PARAMETERS  ARE  USER  DIRECTED 


Display  9  -  Example  Application  of  PLOT  to  a  B-File  Parameter  Time  Histories 
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Display  9  -  Continued 
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5  FREQUENCY  ANALYSIS  PROGRAM  GROUP 


Some  characteristics  of  the  mechanisms  underlying  dynamic  data 
are  not  easily  determined  in  the  time  domain.  It  can  be  more 
appropriate  to  transform  the  parameter  time  history  into  the 
frequency  domain  and  analyze  the  data.  The  frequency  analysis 
program  group  contains  software  which  transforms  time  domain  data 
into  the  frequency  domain.  This  transformation  provides  a  range 
of  analysis  capabilities.  Two  such  capabilities  are  the  determi¬ 
nation  of  the  variation  of  amplitude,  phase,  and  power  spectral 
density  of  a  parameter  time  history  with  frequency  and  the  deter¬ 
mination  of  the  transfer  function  between  two  parameter  time 
histories . 

5.1  Fast  Fourier  Transform  (FFT) 

FFT  performs  a  transformation  of  a  parameter  time  history  into 
the  frequency  domain  through  the  usage  of  the  fast  Fourier 
transform.  The  variation  of  amplitude,  phase,  and  power  spectral 
density  (PSD)  with  frequency  of  a  time  history  is  determined  by 
a  fast  Fourier  transform  of  the  parameter  time  history.  FFT  is 
based  upon  an  algorithm  documented  in  Reference  7  and  implemented 
in  Reference  6. 

5.1.1  Background: 

FFT  performs  a  fast  Fourier  transform  on  a  segment  of  a  parameter 
time  history.  The  fast  Fourier  transform  results  in  a  series  of 
complex  values  known  as  the  spectral  components  of  the  time 
history.  These  spectral  components  occur  at  discrete  frequencies 
and  are  evaluated  at  n/2  +  1  points,  where  n  is  the  length  of  the 
data  sequence.  The  spectral  components  are  complex  with  an 
amplitude  and  phase  (real  and  imaginary  part)  associated  with 
each  spectral  component.  The  amplitude  squared  spectrum  is 
utilized  to  derive  the  power  spectral  density  of  the  sequence 
which  has  the  units  of  energy  per  unit  frequency.  The  amplitude, 
phase,  and  PSD  of  a  parameter  time  history  are  calculated  over  an 
appropriate  frequency  spectrum  which  is  dependent  upon  the  length 
of  the  data  sequence. 

The  fast  Fourier  transform  algorithm  requires  that  the  data 
sequence  to  be  analyzed  consist  of  a  number  of  points  which  is 
equal  to  a  power  of  two.  If  the  input  data  sequence  does  not  meet 
this  requirement  then  two  options  are  available.  The  first  option 
is  to  truncate  the  parameter  time  history  to  a  number  of  points 
equal  to  a  power  of  two  before  analyzing  the  data  sequence.  The 
second  option  is  to  input  the  parameter  time  history  unchanged 
and  FFT  will  zero  fill  the  time  history  to  generate  a  data 
sqeuence  with  a  number  of  points  equal  to  two.  FFT  will 
automatically  check  for  the  power  of  two  criterion  and  proceed 
will  the  zero  fill  if  necessary. 

Leakage  of  spectral  estimates  into  adjacent  frequencies  results 
when  a  time  history  sequence  that  has  been  truncated  is  analyzed 
in  the  frequency  domain.  Data  windows  can  be  employed  to  reduce 


this  leakage.  Most  time  history  data  input  into  FFT  will  be 
truncated  due  to  the  power  of  two  criterion.  To  compensate  for 
this,  FFT  provides  the  option  to  apply  a  Hann  window  to  the  input 
data  sequence  to  reduce  the  impact  of  truncation  effects. 

5.1.2  Input: 

The  parameter  time  history  to  be  analyzed  by  FFT  must  be  resident 
on  a  local  B-file.  This  B-file  is  employed  to  transfer  the  data 
sequence  to  the  analysis  program.  Specifications  required  by  the 
interactive  dialogue  of  FFT  are: 

-  model  number  of  the  Tektronix  graphic  terminal  (if 
employed) 

-  name  of  the  input  B-file 

-  name  of  the  parameter  to  be  analyzed 

-  units  of  the  parameter 

-  data  sample  rate  {samples  per  second] 

-  analysis  start  time-of-day  [seconds] 

-  number  of  points  to  be  analyzed  (power  of  two) 

-  option  to  apply  a  Hann  window 

-  option  to  plot  and/or  list  results 

Depending  on  the  options  selected,  additional  required 
specifications  may  include: 

-  name  of  a  output  file  to  contain  analysis  results 

-  selection  between  graphs  of  amplitude,  phase,  or  PSD  with 
frequency 

-  maximum  value  of  the  plot  ordinate 

-  minimum  value  of  the  plot  ordinate 

Plots  produced  by  FFT  follow  the  conventions  and  format 
associated  with  the  time  history  plot  program  PLOT  discussed  in 
Section  4.2. 


5.1.3  Execution: 

Interactive  execution 
name  call  load  command: 


is  initiated  by  entering  the 


The  program  will  respond  with  a  statement  of  its  capabilities  and 
a  prompt  for  an  interactive  dialogue.  As  indicated  above,  FFT  can 
be  executed  at  a  Tektronix  graphic  terminal  or  a  nongraphic  type 
of  terminal.  Plots  of  the  analysis  results  can  only  be  obtained 
from  a  Tektronix  graphic  terminal.  Beyond  this  restriction 
execution  of  FFT  is  not  dependent  on  the  type  of  terminal 
employed  by  the  user. 


5.1.4  Output: 


FFT  can  be  directed  to  output  analysis  results  to  the  terminal 
screen  and/or  to  a  specified  local  output  file.  FFT  can  also 
display  plots  of  the  analysis  results  on  the  screen  of  a  Tektron¬ 
ix  graphic  terminal.  A  permanent  record  of  the  plotted  results 
can  be  obtained  through  the  usage  of  a  Tektronix  hard  copy  unit. 


5.1.5  Example: 

FFT  was  employed  to  analyze  the  pre-stall  engine  data  in  the 
frequency  domain.  Figure  5  indicates  the  example  application  of 
FFT  to  the  time  history  data  contained  in  B-file  BF024.  Program 
output  consisted  of  a  listing  and  graphic  display  of  the  analysis 
results  on  the  Tektronix  terminal  screen.  Also,  a  listing  of  the 
results  was  directed  to  the  output  file  OF001. 

Display  10  details  the  example  application  of  FFT.  Parameter 
PP128  was  input  through  B-file  BF011.  The  sample  rate  of  PP128 
was  specified  as  1,000  samples  per  second.  The  fast  Fourier 
transform  analysis  was  initiated  at  the  time-of-day  equal  to 
75017.800  seconds  and  included  64  data  points.  This  time  interval 
captured  the  signal  dynamics  of  interest.  The  data  sequence  did 
not  require  zero  filling;  however,  a  Hann  window  was  employed  to 
reduce  leakage  of  the  spectral  estimates  into  adjacent 
frequencies.  A  listing  of  the  analysis  results  was  directed  to 
the  terminal  screen  and  the  output  file.  Display  10  presents  the 
plots  of  the  analysis  results  that  were  displayed  on  the  terminal 
screen. 

The  plots  of  Display  10  show  that  pre-stall  pressure  oscillations 
were  dominant  at  a  frequency  of  60  Hertz.  It  was  suspected  that 
the  resultant  fan  stall  was  caused  by  an  augmentor  combustion 
instability  which  back-pressured  the  fan.  A  combustion  insta¬ 
bility  of  this  type  is  typically  characterized  by  longitudinal 
mode  frequencies  of  50  to  100  Hertz.  Augmentor  rumble  was  further 
suspected  since  the  amplitude  of  the  pressure  fluctuations  were 
measured  to  be  greatest  at  station  5  in  the  fan  duct.  This  is 
shown  by  Figure  9  which  illustrates  a  comparison  of  power 
spectral  density  plots  obtained  at  various  engine  locations. 

A  path  of  propagation  of  the  pressure  oscillations  was  suggested 
by  a  comparison  of  the  amplitudes  of  the  PSD  at  the  characteris¬ 
tic  frequency  of  the  combustion  instability  for  the  various 
measured  flowpath  parameters.  Figure  9  indicates  that  the  source 
of  the  pressure  oscillations  was  in  the  augmentor.  The  pressure 
oscillations  propagated  upstream  in  the  fan  duct.  Upon  reaching 
the  fan  face,  the  pressure  oscillations  propagated  further 
upstream  and  downstream  through  the  engine  core  components. 
However,  the  pressure  oscillations  were  absorbed  by  the  high 
pressure  compressor.  This  resulted  in  stable  pressure  levels  at 
the  core  burner  inlet. 


PT5C  PT2.5C  PT2.5H  PB 
PARAMETER 


Figure  9  -  Power  Spectral  Density  of  Fan  and  Core  Flowpath 
Parameters  Preceding  Engine  Stall 


The  utility  of  the  FFT  analysis  program  can  be  realized  as  a 
primer  for  the  application  of  a  digital  filter  to  enhance  the 
signal  dynamics  of  interest.  FFT  was  employed  for  the  purpose  of 
identifying  the  range  of  frequencies  of  interest  in  the  engine 
pre-stall  dynamics.  As  the  plots  of  Display  10  indicate,  this 
range  of  frequencies  is  from  20  Hertz  to  100  Hertz.  This  informa¬ 
tion  was  utilized  for  the  design  of  the  lowpass  and  multiple  pass 
digital  filters  in  Sections  6.1  and  6.2. 

5.1.6  Limitations  and  nuances: 

FFT  can  accept  a  maximum  of  512  points  from  a  parameter  time 
history. 
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Display  10  -  Example  Application  of  the  Fast  Fourier  Transform  Program  FFT 


RESULTS  OF  FAST  FOURIER  TRANSFORM 
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DO  TOO  WANT  TO  PLOT  THE  FFT  RESULTS’  IY/N) 


Display  10  -  Continued 


POWER  SPECTRAL 


Display  10  -  Continued 


5.2  Transfer  Function  (TRAMS) 

TRANS  computes  the  magnitude  and  phase  of  the  transfer  function 
between  two  parameter  time  histories.  This  capability  is 
particularly  useful  for  determining  relationships  between 
parameters  in  the  frequency  and  time  domain.  TRANS  can  also  be 
used  to  determine  the  frequency  response  of  a  filter.  Such  a 
filter  could  result  from  the  a  design  program  documented  in 
Chapter  6. 

5.2.1  Background: 

TRANS  uses  the  fast  Fourier  transform  algorithm  to  compute  the 
frequency  spectrum  of  the  parameter  time  histories.  Let  the 
frequency  spectrum  of  a  time  history  represented  by  x(n),  where  n 
is  a  point  of  the  data  sequence,  be  given  by  X(w).  Here  w  is  a 
discrete  frequency  in  the  spectrum.  Similarly,  let  the  frequency 
spectrum  of  another  data  sequence  y(n)  be  given  by  Y(w).  Then  the 
transfer  function  between  these  two  parameters  is  represented  as: 

T  (w)  =  Y  ( w )  /X  ( w ) 

T  (w)  is  a  complex  number  whose  value  is  representative  of  the 
transfer  function  between  the  two  parameters  at  all  frequencies 
in  the  spectrum.  Since  T(w)  is  a  ratio  of  spectral  components  and 
can  be  expressed  as  a  vector  with  magnitude  and  phase,  the 
magnitude  of  T(w)  is  calculated  as  201ogT(w)  and  expressed  in 
decibels.  The  phase  is  expressed  in  degrees. 


5.2.2  Input: 

The  input  data  sequences  should  have  a  length  equal  to  a  power  of 
two.  This  restriction  is  applied  because  the  fast  Fourier 
transform  algorithm  is  used  to  compute  the  transfer  function.  If 
the  input  data  sequences  are  not  specified  to  be  analyzed  over  a 
number  of  points  equal  to  a  power  of  two  or  the  data  sequence  is 
not  equal  to  a  power  of  two,  then  the  data  sequence  will  be  zero 
filled  to  the  next  power  of  two. 

A  B-file  format  is  used  to  input  the  two  parameter  time  histories 
for  analysis.  Both  parameter  time  histories  must  be  on  the  B-file 
and  consist  of  the  same  number  of  points.  Specifications  required 
by  the  program  through  the  interactive  dialogue  include: 

-  model  number  of  the  Tektronix  terminal  (if  appropriate) 

-  name  of  the  input  B-file 

-  analysis  start  time-of-day  [seconds] 

-  number  of  data  points  to  be  analyzed 

-  data  sample  rate  [samples  per  second] 

-  name  of  parameter  time  history  1 

-  name  of  parameter  time  history  2 

-  option  to  plot  and/or  list  results 


Depending  on  the  options  selected,  additional  required 
specifications  may  include: 

-  name  of  the  output  file  to  contain  analysis  results 

-  selection  between  a  graph  of  the  amplitude  or  phase  of  the 
transfer  function  versus  frequency 

-  maximum  and  minimum  values  of  plot  ordinate 

Plots  produced  by  TRANS  follow  the  conventions  and  format 
associated  with  the  time  history  plot  program  PLOT  which  is 
discussed  in  Section  4.2. 


5.2.3  Execution: 

Interactive  execution  of  TRANS  is  initiated  by  entering  the 
name  call  load  command: 

TRANS 

The  program  will  respond  with  a  statement  of  its  capabilities  and 
a  prompt  for  interactive  execution.  As  with  the  program  FFT, 
TRANS  can  be  executed  at  a  Tektronix  graphic  terminal  or  a 
nongraphic  type  of  terminal.  Plots  of  the  transfer  function  can 
only  be  obtained  on  a  Tektronix  graphic  terminal.  Besides  this 
restriction,  execution  of  TRANS  is  not  dependent  on  the  type  of 
terminal  employed  by  a  user. 


5.2.4  Output: 

TRANS  can  be  directed  to  list  analysis  results  to  the  terminal 
screen  and/or  to  a  named  local  output  file.  The  listed  analysis 
results  consists  of  a  tabulation  of  the  magnitude  and  phase  of 
the  transfer  function  with  frequency.  This  tabulation  contains 
(n/2)+l  frequency  samples  at  a  spacing  of  fs/n  Hertz,  where  n  is 
the  sequence  length  and  fg  is  the  sample  frequency  in  samples  per 
second. 

A  plot  of  the  transfer  function  magnitude  and  phase  versus 
frequency  can  be  directed  to  the  screen  of  a  Tektronix  graphic 
terminal.  A  permanent  record  of  the  plotted  results  can  be 
obtained  by  a  Tektronix  hard  copy  unit. 


5.2.5  Example: 

Additional  information  concerning  the  mechanisms  which  drive 
dynamic  systems  may  be  obtained  in  the  frequency  domain  by 
determining  the  transfer  function  between  parameter  time 
histories.  As  Figure  5  indicates,  TRANS  was  employed  to  determine 
the  magnitude  and  phase  relationships  of  the  spectral  components 
of  parameters  contained  on  the  input  B-file  BF034.  The  parameters 
contained  on  this  input  file  were  previously  filtered  to  enhance 
the  signal  dynamics.  TRANS  produced  a  listing  and  graphic  display 


of  the  analysis  results  on  the  Tektronix  terminal  screen  as  well 
as  a  listing  of  results  which  was  directed  to  the  output  file 
named  OF002. 

Display  11  details  an  example  application  of  TRANS  to  parameters 
SP128M  and  SP227M.  The  transfer  function  analysis  was  specified 
to  begin  at  a  time-of-day  equal  to  75017.400  seconds  and  to 
encompass  512  data  points.  The  sample  rate  of  the  parameters  was 
1000  samples  per  second.  Parameter  SP128M  was  considered  as  the 
input  signal  and  parameter  SP227M  was  considered  as  the  output 
signal.  A  listing  of  the  analysis  results  was  directed  to  output 
file  OF 0  0 2.  Plots  of  the  transfer  function  magnitude  and  phase 
versus  frequency  were  obtained  at  the  terminal  screen  and  are 
given  in  Display  11. 

TRANS  is  particularly  useful  in  determining  the  frequency 
response  of  a  digital  filter.  In  this  case,  the  two  signals 
analyzed  are  the  input  and  output  of  the  digital  filter.  An 
example  of  this  application  of  TRANS  is  presented  in  Section  7.1. 

5.2.6  Limitations  and  nuances: 

TRANS  can  accept  a  maximum  of  512  points  from  a  parameter  time 
history. 


THIS  PP0GRAI1  COMPUTES  A  TRANSFER  function  between  two  signal 
SEQUENCES,  the  input  B-FlLE  HOST  CONTAIN  ThC  INPUT  AND  OUTPUT 
SEQUENCES.  OUTPUT  IS  USER  DIRECTED  TO  BE  LISTED  OR  PLOTTED 
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Display  11  -  Example  Application  of  the  Transfer  Function  Propram  TRANS 
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6.  DIGITAL  FILTER  PROGRAM  GROUP 

The  digital  filter  program  group  provides  the  capability  to 
design,  analyze,  and  implement  digital  filters  on  a  parameter 
time  history.  The  software  comprising  this  program  group  are 
modified  versions  of  finite  impulse  response  (FIR)  filter  design 
and  implementation  programs.  These  programs  were  obtained  from 
the  digital  signal  processing  software  package  which  is 
documented  in  Reference  2.  The  digital  signal  processing  software 
package  provides  additional  filtering  capability;  however,  the 
filter  design  and  implementation  programs  discussed  below  will 
meet  most  filtering  requirements. 

Much  of  the  background  discussion  presented  below  and  in  the 
following  sections  is  paraphrased  from  Reference  2.  A  more  de¬ 
tailed  discussion  of  digital  filters  is  given  in  this  reference. 
Additionally,  an  overview  of  digital  signal  processing 
terminology  is  given  in  Reference  8. 

The  frequency  spectrum  that  a  signal  encompasses  often  represents 
the  characteristics  of  the  signal  more  explicitly  than  its  time 
domain  counterpart.  When  analyzing  a  signal  in  the  frequency  do¬ 
main,  spectral  components  are  revealed  that  may  be  removed  or 
enhanced  by  applying  a  digital  filter.  An  example  of  this  situa¬ 
tion  is  noise  due  to  a  power  supply  which  is  at  a  discrete 
frequency.  It  would  be  desirable  to  eliminate  this  noise  from  the 
otherwise  good  signal.  This  can  be  accomplished  by  suppressing 
the  components  of  the  frequency  spectrum  which  correspond  to  the 
frequency  of  the  power  supply  noise.  This  capability  is  provided 
by  the  filter  design  programs  discussed  in  this  chapter. 

The  design  of  a  digital  filter  is  determined  by  a  number  of 
specifications  which  describe  the  intended  function  of  the 
filter.  These  specifications  designate  the  various  pass  bands  and 
stop  bands  of  the  filter.  The  band  of  frequencies  that  are  to 
pass  essentially  unattenuated  through  the  filter  is  called  the 
pass  band.  Similarly,  the  band  of  frequencies  that  are  to  be 
attenuated  is  called  the  stop  band.  A  transition  band  is  a  range 
of  frequencies  which  the  filter  will  partially  attenuate  the 
frequency  components.  Although  it  is  desirable  to  not  specify  a 
transition  band,  this  is  not  obtainable  for  a  real  filter. 
Various  combinations  of  the  stop,  pass,  and  transition  bands 
produce  typical  filter  designs  as  illustrated  in  Figure  10. 

Within  each  pass  band  and  stop  band,  the  affect  of  a  filter  may 
be  characterized  by  the  amplitude  and  phase  response  of  the 
filter.  The  amplitude  response  of  the  filter  is  the  ratio  of  the 
output  signal  to  the  input  signal.  The  amplitude  response  is 
expressed  in  decibels  at  any  component  frequency.  The  phase 
response  of  the  filter  indicates  the  time  delay  of  the  output 
signal  relative  to  the  input  signal.  Phase  response  is  expressed 
in  degrees  at  any  component  frequency.  The  amplitude  and  phase 
response  specify  the  frequency  response  of  the  filter. 


Thefilter  design  programs  contained  in  the  ISAS  produce  a 
finite  impulse  response  (FIR)  filter.  An  FIR  filter  produces  a 
linear  phase  response  with  frequency.  This  corresponds  to 
imparting  a  constant  time  delay  to  the  filter  output  time  history 
relative  to  the  input  time  history.  However,  this  character istic 
can  be  compensated  through  a  temporal  alignment  of  the  parameter 
time  history.  This  can  be  accomplished  by  using  the  data 
conditioning  program  TIMALIN. 

The  above  discussion  provides  a  basic  understanding  of  digital 
filters.  This  information  and  that  contained  in  the  following 
sections  should  facilitate  implementation  of  a  digital  filter 
with  a  reasonable  level  of  confidence.  Again,  a  more  detailed 
discussion  of  digital  filters  is  given  in  References  2  and  8. 


6.1  Lowpass  Filter  Design  (LOWPASS) 

The  filter  design  program  LOWPASS  generates  a  filter  file  for  a 
maximally  flat  pass  band  and  stop  band  FIR  lowpass  filter.  The 
resultant  filter  file  is  used  by  the  filter  implementation 
program  to  filter  a  parameter  time  history. 


6.1.1  Background: 

LOWPASS  is  characterized  by  an  ideal  magnitude  response  that  has 
a  value  of  1  within  the  pass  band  and  a  value  of  0  within  the 
stop  band.  Since  the  pass  band  and  stop  band  cannot  overlap,  a 
transition  band  is  employed.  This  transition  band  determines  the 
steepness  of  the  filter  rolloff  and  the  order  of  the  filter. 
LOWPASS  solves  for  a  closed  form  expression  of  the  filter 
response  which  yields  filter  coefficients.  These  filter 
coefficients  comprise  the  filter  description  that  is  written  to 
the  filter  file. 


6.1.2  Input: 

Input  to  the  filter  design  program  LOWPASS  is  accomplished  only 
through  an  interactive  dialogue  with  the  program.  Specifications 
required  for  the  filter  design  are: 

-  sample  rate  of  the  data  to  be  filtered  [samples  per 
second] 

-  center  frequency  of  the  transition  band  (less  than  1/2  the 
sample  frequency)  [Hertz] 

-  transition  band  frequency  width  (greater  than  the  sample 
frequency  /  25)  [Hertz] 

-  option  to  write  the  filter  description  to  a  filter  file 

-  option  to  list  the  filter  coefficients  at  the  terminal 
screen  or  to  an  output  file 

Depending  on  the  options  selected  above,  additional  required 
specifications  may  include: 


-  name  for  the  filter  file 

-  name  for  the  output  file 

Upon  completion  of  one  filter  design,  another  lowpass  filter  may 
be  designed  by  reinitiating  the  program  inquiry  and  providing  the 
above  specifications  again. 

6.1.3  Execution: 

Interactive  execution  of  LOWPASS  is  initiated  by  entering  the 
name  call  load  command: 

LOWPASS 

The  program  will  respond  with  a  statement  of  its  capabilities  and 
a  prompt  for  an  interactive  dialogue.  Upon  termination  of 
LOWPASS,  the  terminal  is  returned  to  the  command  level  of 
NOS. 


6.1.4  Output: 

LOWPASS  produces  a  filter  description  file  which  contains  the 
filter  coefficients.  This  file  is  input  to  the  filter 
implementation  program  FILTER.  Also,  the  filter  description  can 
be  listed  at  the  terminal  screen  or  written  to  a  local  output 
file  for  saving  on  a  permanent  storage  device  or  listing  at  a 
printer . 

6.1.5  Example: 

Digital  data  filtering  techniques  were  applied  to  the  parameter 
time  histories  obtained  during  transient  augmentor  operation. 
The  filter  pass  band  was  specified  to  pass  a  band  of  frequencies 
which  contained  the  characteristic  frequency  of  the  pre-stall 
flowpath  pressure  dynamics.  Lowpass  provided  a  filter  design 
which  would  enhance  the  pressure  oscillation  data.  Figure  5 
indicates  the  application  of  LOWPASS  to  generate  the  filter  file 
FF001  which  was  subsequently  used  by  FILTER.  Also,  a  listing  of 
the  filter  description  was  directed  to  the  local  output  file 
OF003. 

Display  12  provides  a  record  of  an  interactive  filter  design 
session  which  exemplifies  the  dialogue  and  execution  of  LOWPASS. 
The  sample  rate  of  the  data  to  be  filtered  was  specified  as  1,000 
samples  per  second.  From  the  results  of  the  fast  Fourier 
transform  analysis,  it  was  determined  that  the  characteristic 
frequency  of  the  engine  pre-stall  dynamics  was  below  200  Hertz. 
Thus,  the  center  frequency  of  the  transition  band  width  was 
specified  as  250  Hertz  and  the  transition  band  width  was 
specified  to  be  60  Hertz.  The  resultant  filter  description  was 
written  to  the  filter  file  named  FF001.  Display  12  also  shows  the 


screen  listing  of  the  filter  coeffifients  that  were  directed  to 
the  output  file  OF003. 

Application  of  the  filter  design  to  the  parameter  time  histories 
is  discussed  in  Section  6.3.  Also,  the  LOWPASS  filter  design  was 
verified  by  utilizing  the  white  noise  data  file  to  drive  the 
input  of  the  filter  design.  The  output  of  the  filter  thus 
indicated  the  frequency  response  of  the  filter.  This  procedure  is 
presented  in  Chapter  7. 

6.1.6  Limitations  and  nuances: 

Some  restrictions  are  imposed  upon  the  design  specifications  for 
LOWPASS.  These  restrictions  are: 

1.  The  center  frequency  of  the  transition  band  must  be 
greater  than  the  sample  frequency  divided  by  2. 

2.  The  width  of  the  transition  band  must  be  greater  than  the 
sample  frequency  divided  by  25. 

If  these  restrictions  are  observed,  LOWPASS  will  normally 
converge  on  a  filter  design. 


Display  12  -  Example  Utilization  of  LOWPASS  for  a  Digital  Filter  Design 
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6.2  Multiple  Passband  Filter  Design  (MLTPASS) 

MLTPASS  designs  a  finite  impulse  response  filter  based  upon  given 
specifications . 


6.2.1  Background: 

Many  filter  configurations  may  be  designed  by  using  the  multiple 
pass  band  and  stop  band  specifications  of  MLTPASS.  The  amplitude 
response  of  any  filter  is  characterized  by  at  least  one  pass  band 
and  one  stop  band  with  a  transition  band  between  these  two  bands. 
The  lowpass  and  highpass  filters  are  characterized  by  one  pass 
band  and  one  stop  band.  However,  other  filter  designs  may  have 
more  than  one  stop  band  and  one  pass  band.  For  example,  bandpass 
and  band  reject  filters  have  two  pass  bands  and  stop  bands. 
Figure  10  indicates  the  relationship  between  the  pass  band  and 
the  stop  band  frequencies  for  various  types  of  filters.  The  need 
for  the  specification  of  the  lower  and  upper  cutoff  frequencies 
as  well  as  the  the  lower  and  upper  stop  band  edge  frequencies  is 
dependent  upon  the  type  of  filter. 

MLTPASS  uses  an  Remez  exchange  algorithm  to  determine  filter 
values  which  minimize  a  weighted  Chebyshev  error  function  over  a 
grid  of  frequencies.  This  means  that  the  difference  between  the 
desired  filter  function  and  a  trial  filter  function  is  calculated 
over  a  set  of  frequencies.  The  magnitude  of  this  difference  is 
then  used  to  update  the  filter  function.  Also,  the  extent  to 
which  error  influences  the  filter  function  in  each  of  the 
specified  bands  is  determined  by  the  weighting  value  for  that 
band.  Thus,  the  value  of  the  band  weighting  signifies  the 
relative  importance  of  the  filter  response. 


6.2.2  Input: 

Input  to  MLTPASS  is  provided  through  interactive  dialogue  with 
the  program.  Specifications  required  by  MLTPASS  are: 

-  filter  length,  i.e.  the  number  of  points  to  comprise  the 
FIR  filter  (maximum  is  128) 

-  number  of  frequency  bands  (disjoint  pass  bands  and 
stop  bands) 

-  grid  density  (nominal  =  16) 

-  sample  rate  of  the  data  to  be  filtered  [samples  per 
second ] 

-  edge  frequencies  of  the  passbands  and  stop  bands  [Hertz] 

-  filter  gain  (1  for  an  ideal  passband  or  0  for  an  ideal 
stopband ) 

-  relative  weighting  for  each  specified  band 

-  option  to  write  the  filter  description  to  a  filter  file 

-  option  to  list  the  filter  coefficients  at  the  terminal 
screen  or  to  an  output  file 

-  gain  of  the  filter  for  each  specified  frequency  band 
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Depending  on  the  options  selected  above,  additional  required 
specifications  may  include: 


-  name  for  the  filter  file 

-  name  for  the  output  file 

Upon  completion  of  one  filter  design,  another  filter  may  be 
designedby  responding  affirmatively  to  reinitiating  the  program 
inquiry  and  providing  the  appropriate  specifications  from  those 
given  above. 


6.2.3  Execution: 

Interactive  execution  of  MLTPASS  is  initiated  by  entering  the 
name  call  load  command: 

MLTPASS 

The  program  will  respond  with  a  statement  of  its  capabilities  and 
a  prompt  for  an  interactive  dialogue.  The  terminal  is  returned  to 
the  command  level  of  NOS  upon  termination  of  program  execution. 


6.2.4  Output: 

MLTPASS  produces  a  filter  description  file  which  contains  the 
filter  coefficients.  This  file  is  input  to  the  filter 
implementation  program  FILTER  to  actually  filter  a  parameter  time 
history.  This  filter  description  file  can  be  listed  at  the 
terminal  screen  or  written  to  a  local  output  file. 


6.2.5  Example: 

The  MLTPASS  filter  design  program  was  utilized  to  design  a 
digital  filter  which  would  enhance  the  pre-stall  dynamic  engine 
data.  Figure  5  illustrates  the  utilization  of  MLTPASS  for  the 
generation  of  the  filter  file  FF002  which  was  subsequently  used 
by  FILTER  to  filter  the  parameter  time  histories  of  interest. 

Display  13  provides  a  record  of  an  interactive  filter  design 
session  which  exemplifies  the  execution  and  dialogue  of  MLTPASS. 
A  multiple  pass  band  and  stop  band  FIR  digital  filter  was 
designed  for  the  enhancement  of  data  obtained  during  dynamic 
oscillations  of  the  gas  turbine  engine.  The  data  to  be  filtered 
is  that  which  exemplified  the  execution  of  the  Fourier  analysis 
program  FFT  and  the  usage  of  the  filter  design  program  LOWPASS. 
Analysis  of  the  dynamic  data  using  FFT  further  determined  the 
characteristic  frequency  of  the  pre-stall  dynamics  to  be  centered 
about  a  frequency  of  60  Hertz.  The  MLTPASS  filter  design  was 
specified  to  generate  a  pass  band/stop  band  filter  to  further 
enhance  the  pre-stall  dynamics.  The  pass  band  was  specified  to 
range  from  0  Hertz  to  100  Hertz  and  the  stop  band  was  specified 
to  range  from  120  Hertz  to  500  Hertz.  500  Hertz  was  chosen  as  the 


upper  cutoff  frequency  because  this  is  the  folding  frequency 
corresponding  to  the  data  sample  rate  of  1000  samples  per  second. 
This  filter  design  resulted  in  a  transition  band  from  100  Hertz 
to  120  Hertz.  The  resultant  filter  design  was  saved  on  the  filter 
file  named  FF002  for  subsequent  application  through  the  filter 
implementation  program  FILTER. 

Additional  specifications  were  input  to  MLTPASS.  The  number  of 
filter  coefficients  was  chosen  to  be  64.  Since  one  stop  band  and 
one  pass  band  were  desired,  the  total  number  of  frequency  bands 
specified  was  two.  The  filter  grid  density  was  set  at  the  nominal 
value  of  16.  The  sample  rate  of  the  parameter  time  history  to  be 
filtered  was  given  as  1,000  samples  per  second.  The  gain  of  the 
pass  band  and  the  stop  band  were  set  equal  to  1  and  0 
respectively.  The  frequency  bands  were  given  a  relative  weighting 
of  10.  A  listing  of  the  resultant  filter  description  was  directed 
to  the  terminal  screen.  Also,  output  was  directed  to  an  output 
file  named  OF004  for  subsequent  listing  at  a  line  printer. 

Verification  of  the  filter  design  produced  by  MLTPASS  was 
accomplished  by  utlizing  the  white  noise  data  file.  This  file  was 
used  as  input  to  the  filter  design  produced  above.  The  resultant 
filter  output  indicated  the  frequency  response  of  the  filter. 
This  verification  process  is  presented  in  Chapter  7. 


6.2.6  Limitations  and  nuances: 


each 
e  of 


importance  of  the  filter  response. 


Display  13  -  Example  Utilization  of  MLTPASS  for  a  Digital  Filter  Design 
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6.3  Filter  Implementation  (FILTER) 

FILTER  filters  a  parameter  time  history  by  using  a  previously 
generated  filter  design. 

6.3.1  Background: 

The  filtering  of  a  parameter  time  history  is  accomplished  by  the 
convolution  of  the  filter  design  with  the  data  sequence.  This 
convolution  is  accomplished  by  using  a  fast  Fourier  transform. 
The  fast  Fourier  transform  algorithm  requires  that  the  data 
sequence  to  be  filtered  consist  of  a  number  of  points  which  is 
equal  to  a  power  of  two.  If  the  input  data  sequence  does  not  meet 
this  requirement  then  FILTER  will  zero  fill  the  input  data 
sequence  to  the  next  power  of  two. 

6.3.2  Input: 

A  B-file  format  is  employed  to  input  the  parameter  time  history 
to  FILTER.  Also,  a  filter  file  which  has  been  previously 
generated  is  used  to  relay  the  particular  filter  design 
specifications  to  the  program.  Specifications  requested  by  FILTER 
through  an  interactive  dialogue  are: 

-  nameof  the  B-file  which  contains  the  data  sequence  to 
be  filtered 

-  name  for  the  B-file  generated  by  FILTER  which  will  contain 
the  filtered  parameter  time  history 

-  name  of  the  parameter  to  be  filtered 

-  name  for  the  filtered  parameter 

-  selection  of  terminal  entry  or  a  filter  file  as  thesource 
of  the  filter  description 

If  a  terminal  entry  was  selected  then  additional  requested 
specifications  are: 

-  length  of  the  filter 

-  filter  sequence  point  number  and  gain  for  each  non-zero 
filter  value 

If  a  filter  file  was  selected  as  the  filter  description  source 
then  an  additional  requested  specification  is: 

-  name  of  the  filter  file 

If  it  is  desired  to  filter  another  parameter  time  history,  the 
above  specifications  must  be  entered  again. 

6.3.3  Execution: 

Interactive  execution  of  FILTER  is  initiated  by  entering  the 
name  call  load  command: 


FILTER 


The  program  will  respond  with  a  statement  of  its  capabilities  and 
a  prompt  for  interactive  execution.  Termination  of  FILTER  returns 
the  terminal  to  the  command  level  of  NOS. 


6.3.4  Output: 

FILTER  produces  a  filtered  parameter  time  history  which  is 
appended  to  the  input  time  history  and  output  on  a  B-file. 


6.3.5  Ezaaple: 

* 

i  FILTER  was  employed  to  enhance  the  dynamic  parameters  of  interest 

through  filter  designs  supplied  by  the  design  programs  LOWPASS 
and  MLTPASS.  Figure  5  illustrates  the  application  of  FILTER  to 
the  parameter  time  histories  contained  on  B-file  BF011.  In  this 
example,  FILTER  was  employed  twice,  once  with  the  LOWPASS  filter 
design  given  by  FF001  and  again  with  the  MLTPASS  filter  design 
given  by  FF002.  The  resultant  filtered  parameters  were  output  on 
B-file  BF024. 

Display  14  provides  an  example  application  of  the  LOWPASS  filter 
design  to  the  time  history  of  PP128.  The  filtered  parameter  time 
history  FP128L  was  directed  to  an  intermediate  B-file  named 
•  BF021.  The  result  of  filtering  PP128  is  illustrated  in  Display  14 

1  by  comparing  the  plots  of  the  PP128  and  FP128L.  Qualitatively, 

the  filtered  parameter  shows  a  reduction  of  the  higher  frequency 
v  oscillations  which  reflects  the  stopband  specifications  of  the 

)  LOWPASS  filter  design.  Also,  the  LOWPASS  filter  design  does  not 

impart  a  time  delay  to  the  filtered  time  history.  Parameter 
j  FP128L  is  temporally  aligned  with  the  input  time  history  PP128. 

FILTER  was  also  applied  to  parameter  PP227  through  the  multiple 
bandpass  filter  design  given  by  MLTPASS.  The  result  of  this 
filter  application  is  shown  in  Display  14  by  the  plot  of  the 
filtered  parameter  FP227M.  This  plot  can  be  compared  to  the  plot 
of  the  input  time  history  of  parameter  PP227.  As  the  plot  of 
FP227M  indicates.  The  MLTPASS  filter  design  imparts  a  time  delay 
to  the  filtered  parameter  time  history.  This  time  delay  may  be 
corrected  through  the  application  of  TIMALIN  as  given  in  Section 
3.5. 

The  filtering  of  a  time  history  by  FILTER  is  accomplished  by  a 
convolution  of  the  filter  design  with  the  data  sequence.  Since 
the  filter  has  a  finite  length  N,  N/2  end  points  on  each  side  of 
the  filtered  time  history  will  be  inaccurate.  This  is 
dramatically  shown  by  the  plot  of  FP227M.  In  this  case,  this 
effect  is  combined  with  the  time  lag  due  to  the  linear  response 
of  the  filter  design  by  MLTPASS. 


The  digital  filter  may  uniformly  attenuate  the  level  of  the 
signal.  This  characteristic  is  most  pronounced  by  the 


implementation  of  the  lowpass  filter  as  is  illustrated  by  FP128L. 
A  good  estimation  of  the  expected  attenuation  can  be  determined 
by  filtering  a  time  history  of  white  noise.  A  white  noise  data 
file  is  provided  in  the  fileset  for  this  purpose  and  is  discussed 
in  Section  7.1.  The  filtered  white  noise  can  then  be  compared  to 
the  original  white  noise  in  the  frequency  domain.  This  will 
determine  the  expected  attenuation  of  a  signal  filtered  with  the 
filter  design  over  a  range  of  frequencies  of  interest. 

The  digital  filtering  techniques  described  above  were  applied  to 
the  dynamic  parameters  measured  prior  to  engine  stall.  Figure  11 
illustrates  the  digital  filter  processing  results  for  the 
flowpath  parameters  PS5C,  PT2.5H,  PB,  and  the  augmentor  fuel 
manifold  pressure  PSEG5.  These  parameters  were  filtered  in  a 
frequency  band  from  40  Hertz  to  70  Hertz  and  over  a  period  of 
time  which  encompassed  the  dynamics  of  the  fuel  system  and  the 
core  and  bypass  flowpaths.  This  period  of  time  corresponds  to 
elapsed  times-of-day  from  75017.200  seconds  to  75018.200  seconds. 
The  pressure  oscillations  of  PSEG5  are  clearly  shown  and  are 
observed  to  influence  the  fan  duct  pressure  PS5C  and  the  core 
flowpath  pressure  PT2.5H.  This  influence  was  probably  due  to  a 
combustion  instability  driven  by  oscillating  fuel  delivery  rates 
to  the  augmentor. 

Digital  filtering  techniques  were  employed  to  determine  if  the 
characteristic  frequency  of  the  augmentor  instability  was 
constant  for  the  period  of  time  immediately  preceding  stall. 
Figure  12  shows  the  results  from  filtering  PS5C  into  four 
frequency  bands  which  range  from  30  Hertz  to  70  Hertz.  During  the 
initial  growth  of  the  pressure  oscillation,  the  higher  frequency 
bands  contained  most  of  the  frequency  content.  However,  as  the 
stall  point  approached,  the  characteristic  frequency  of  the 
flowpath  oscillation  was  shifted  to  the  lower  frequency  bands 
where  it  was  most  pronounced  before  stall. 

6.3.6  Limitations  and  nuances: 

The  end  points  of  the  filtered  parameter  time  history  will  be 
inaccurate  due  to  the  finite  length  of  the  filter.  The  number  of 
end  points  affected  is  dependent  upon  the  length  of  the  filter. 


EQ6  PS5C  PT2.5H 

SI)  (PSI)  (PSD 


PS5C  PS5C  PS5C  PS5C 
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Figure  12  -  Multiple  Frequency  Band  Filter  Results 

Prior  to  Engine  Stall 
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Example  Application  of  FILTER  Using  L0VJPASS 
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7.  TEST  SOURCE  DATA  FILES 


Test  source  data  files  provide  a  means  to  test  the  output  of  an 
analysis  program  with  a  well  defined  input.  Currently,  one  such 
data  file  is  provided  on  the  ISAS  fileset. 


7.1  White  Noise  (DATA/WHITE) 

White  noise  is  represented  by  a  data  sequence  on  the  file 
named  DATA/WHITE.  This  file  is  an  UFTAS  standard  B-file  with  the 
fileset  group  name  DATA  and  element  name  of  WHITE.  The  white 
noise  file  provides  a  source  test  signal  with  frequency  content 
spread  across  the  sampled  data  spectrum.  This  signal  may  be  used 
to  drive  the  input  of  a  digital  filter  under  consideration  so 
that  the  filter's  output  is  an  indication  of  the  frequency 
response  of  the  filter. 


7.1.1  Background: 

After  a  filter  design  has  been  generated,  it  is  useful  to  examine 
the  actual  filter  response  and  determine  if  the  actual  filter 
response  meets  the  design  specifications.  A  simple  way  to 
accomplish  this  is  to  provide  a  broad  spectrum  time  history  as 
input  to  the  filter  design  and  use  a  spectrum  analysis  tool  to 
examine  the  output  of  the  filter.  Such  a  broad  spectrum  input  is 
white  noise  which  contains  sinusoids  of  equal  amplitude  over  a 
broad  frequency  range. 

The  fast  Fourier  transform  may  be  used  to  analyze  the  frequency 
response  of  a  system  subjected  to  a  particular  signal.  If  a 
digital  signal  is  composed  of  sinusoids  evenly  spaced  across  the 
frequency  spectrum,  the  FFT  will  calculate  frequency  content 
within  evenly  spaced  intervals  across  the  spectrum.  Also,  if  the 
N-point  FFT  of  an  sequence  is  calculated  from  the  signal,  the 
frequency  content  at(N/2)+l  positive  frequencies  is  obtained  by 
the  fast  Fourier  transform.  These  frequencies  are  located  between 
0  and  fs/2  and  have  a  spacing  of  fs/N,  where  fs  is  the  sample 
frequency. 

If  the  sinusoids  of  the  test  signal  are  located  at  precisely  the 
center  frequencies,  no  smearing  of  the  frequency  spectrum  will 
occur.  The  frequency  content  of  the  signal  will  be  accurately 
represented  by  a  FFT.  Correspondingly,  the  output  of  a  linear 
system  driven  by  this  input  will  also  have  frequency  content  at 
only  the  center  frequencies.  The  output  of  the  system  may  be 
analyzed  to  determine  the  effectiveness  of  the  system  on  the 
input  signal.  Thus  the  filtered  output  from  a  digital  filter  that 
had  as  input  the  white  noise  data  file  may  be  compared  to  this 
input  to  determine  the  effectiveness  of  the  digital  filter.  This 
comparison  is  best  conducted  in  the  frequency  domain  by  employing 
a  frequency  analysis  program  such  as  TRANS  or  FFT. 


7.1.2  Data  file  specifications: 


Specifications  of  the  white  noise  data  file  are  given  below  to 
facilitate  the  application  of  this  test  sequence  in  the  digital 
filter  and  frequency  analysis  programs. 


File  name 
File  type 
Parameter  name 

Data  sample  rate 
Start  time-of-day 
Finish  time-of-day 
Signal  length 
Number  of  sinusoids 
Sinusoid  spacing 


DATA/WHITE 

B-file 

N0ISE1  (this  is  the  only  parameter  on 
the  B-file) 

1,000  samples  per  second 

0.001  seconds 

1.024  seconds 

1024  samples 

513  (including  DC) 

0.97656  Hertz  between  0  and  500  Hertz 


7.1.3  Example: 

The  data  file  named  WHITE  was  employed  to  verify  the  digital 
filter  designs  detailed  in  Sections  6.1  and  6.2.  This  white  noise 
B-file  was  generated  through  the  UFTAS  utility  data  generation 
program  DATAGEN  which  is  documented  in  Reference  5.  Depending  on 
the  specifications  of  the  filter  design  (e.g.  data  sample  rate, 
frequency  band  range,  etc.),  it  may  be  necessary  to  generate  a 
white  noise  data  file  which  reflects  the  sample  rate  and 
frequency  content  of  the  data  to  be  filtered.  The  white  noise 
data  specifications  given  above  were  specified  through  DATAGEN  to 
be  compatible  with  the  example  test  data.  For  reference,  a  plot 
of  WHITE  is  given  by  Display  15  which  was  obtained  from  PLOT. 

As  Figure  5  indicates,  WHITE  was  employed  as  input  to  the  filter 
implementation  program  FILTER  to  verify  the  digital  filter 
designs.  The  parameter  N0ISE1  was  filtered  according  to  the 
LOWPASS  design  specifications  input  through  the  filter  file 
FF 0  01.  The  filtered  parameter  time  history  named  LFNOIS  was 
output  on  the  B-file  WN024.  This  filtered  time  history  was 
subsequently  compared  with  the  unfiltered  time  history  N0ISE1  in 
the  frequency  domain  by  TRANS.  Display  16  shows  the  resultant 
plot  of  the  transfer  function  magnitude  of  the  transfer  function 
between  these  two  parameters.  The  effect  of  the  LOWPASS  filter 
design  is  evident  and  corresponds  to  the  filter  specifications. 
The  pass  band  is  present  from  0  Hertz  to  200  Hertz  with  the 
amplitude  of  the  frequency  components  above  250  Hertz  severely 
attenuated  by  the  filter.  However,  Display  16  shows  that  the  pass 
band  frequency  components  are  slightly  attenuated  by 
approximately  2  dB.  This  has  been  observed  to  be  a  consequence  of 
the  LOWPASS  filter  design  and  should  be  noted  for  reference. 

In  a  similar  manner,  the  white  noise  parameter  NOISEl  was 
employed  to  verify  the  MLTPASS  filter  design.  Display  17  shows  a 
plot  of  the  transfer  function  magnitude  between  NOISEl  and  the 
filtered  white  noise  parameter  SFNOISE.  The  design  specifications 
of  the  filter  file  FF002  are  reflected  in  this  plot  of  the 


transfer  function  magnitude.  The  pass  band  is  predominant  from  0 
Hertz  to  100  Hertz  without  any  attenuation.  The  transition  band 
displays  rapid  amplitude  dropoff  of  the  filtered  parameter  by 
approximately  25  dB.  Also,  the  stop  band  is  severely  attenuated 
to  the  limiting  frequency  of  500  Hertz. 

Verification  of  the  filter  designs  employed  to  enhance  the  pre¬ 
stall  test  data  through  the  above  procedure  provided  confidence 
for  the  implementation  of  the  filter  designs  to  the  test  data. 
Also,  insight  was  gained  into  the  exact  behavior  of  the  filter 
over  the  spectral  components  of  the  white  noise  time  history. 
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ANALYSIS  Of  GAS  TURBINE  ENGINE  DYNAMIC  INSTABILITIES 


Ronald  M.  Devlno* 

Air  Forca  Flight  Taat  Cancer 
Edvarda  Air  Foret  Bata,  California 


Abatract 

A  aachod  hat  baan  developed  for  the  analysis  of 
dynamic  lnacabllltlaa  that  have  baan  obeerved  In 
gaa  turbine  englnaa  during  flight  teat.  An  inves¬ 
tigation  which  ezaapllflea  data  analysis  techniques 
waa  conducted  on  dynamic  data  recorded  during  »ug- 
aentor  cranaiene  operation.  Signal*  from  preaaure 
tranaducara  located  In  the  angina  gaa  path  and 
fuel  delivery  ayataa  ware  recorded  and  digitized. 
Oaca  aaalyala  techniques  Included  ciaa  history 
plotting  of  low  and  high  frequency  raaponaa  data, 
Fourier  Transform  baaed  spectral  analysis,  and 
digital  filtering.  Pre-stall  dynamic  Instabilities 
ware  characterized  by  low  frequency  oscillations 
throughout  many  angina  components.  The  anal¬ 
ysis  of  component  lnceracdons  Identified  the 
path  of  dynamic  Instabilities  through  the  angina. 
A  transfer  of  dynamic  energy  fraa  one  frequency 
band  lnco  another  was  also  Identified. 


Nomenclature 

AC  type  of  oscillating  signal  employed  by 

masking  technique 

AJ  exhaust  noazla  area,  percent 

EPR  engine  pressure  redo 

EU  engineering  units 

FPT  Fast  Fourier  Transform 

EM  frequency  modulated  signal 

LOO  augment or  light-off  detector 

MASK  type  of  constant  slaaal  employed  by  mask¬ 
ing  technique 

Ml  fan  rotor  speed,  tpm 

M2  core  rotor  speed,  rpa 

OKIG  type  of  dynamic  signal  employed  by  mask¬ 
ing  technique 

Pb  core  burner  lnlac  prasaure,  lbf/ln2 

PCM  pulse  code  module  tad  signal 

PLA  power  lever  angle 

Psegl  fuel  preaaure  In  augmentor  segment  1, 
lbf/ln2 


Psegl  fuel  prasaure  In  augmentor  segment  2, 
lbf/in2 

Psegl  fuel  pressure  In  augmentor  segment  3, 
lbf/in2 


•Aerospace  Engineer,  Flight  Test  Technology  Branch, 
Member  ALLA 


PsagA  fuel  prasaure  in  augmentor  secant  4, 
lbf/ln2 

PsagS  fuel  prasaure  in  augsentor  segment  3, 
lbf/ln2 

Ps2.Sc  fan  exit  duct  static  pressure,  lbf/ln2 

Pt2.5c  fan  exit  duct  total  pressure,  lbf/ln2 

Pt2.Sh  fan  exit  core  total  pressure,  lbf/ln2 

PsSc  fan  duct  static  pressure,  lbf/ln2 

Pt6m  turbine  discharge  (mixed  core  and  fan 
scream)  total  prasaure,  lbf/ln2 

SANDS  Standard  Analog  to  Digital  System  for 
processing  data 

c  elapsed  time,  seconds 

UFTAS  Uniform  Flight  Test  Analysis  System 


Introduction 

A  continuous  technological  challenge  exists  for 
high  performance  aircraft  capable  of  executing 
multiple  mission  roles.  This  challenge  has  driven 
the  design  of  the  aircraft  gas  turbine  engine 
throughout  Its  history  of  development.  Today's 
typical  high  performance  gaa  turbine  engine  Is  an 
augmented  turbofan  engine  with  features  such  as 
variable  compresaor  vane  geometry,  variable  exhaust 
nossle  area,  and  a  hydromechanical  or  electronic 
control  system.  Such  a  design  provides  a  high 
chrust-co-welght  ratio  and  optlmia  performance  over 
a  wide  range  of  operating  conditions.  Continuing 
gas  turbine  technological  advancements  coupled  with 
aircraft  Integration  through  the  engine  control 
systms  have  resulted  In  complex  designs  to  meet 
operational  requirements. 

The  complex  design  of  the  modern  gas  curblne 
engine  often  harbors  operational  anomalies. 
Operation  of  the  aircraft’s  propulsion  system  at 
extreme  regions  of  the  flight  envelope  can  cause 
components  of  the  gas  curblne  engine  to  operate  off 
their  design  point.  Components,  such  as  che  high- 
pressure  compressor,  can  experience  rotating  stall 
or  surge  which  adversely  affects  other  engine 
components.  The  transient  processes  which  occur 
during  a  rapid  advance  of  che  chroctle  from  an 
idle  or  military  power  setting  to  a  maxima 
power  setting  may  Induce  combustion  Instabilities 
in  che  afterburner  and  its  control  system.  Combus¬ 
tion  instabilities  often  inherent  in  high  perfor¬ 
mance  thrust  augmentation  systems  can  result  in 
Large  amplitude  pressure  oscillations.  These 
pressure  oeclllacions  can  stall  che  compression 
components  of  che  engine  or  Induce  augmencor  blow¬ 
out.  Transient  processes  resulting  from  a  throttle 
excursion  can  also  impose  unaeecable  demands  on 
che  engine  concrol  systms.  Phenomena  and  processes 


of  Choao  typo*  within  tha  |ta  turbine  angina  re¬ 
strict  oparatlon  of  tha  aircraft  ovar  a  ragloa  of 
its  flight  anvalopa  and  raduca  tha  affaetlvanaaa 
of  lta  alaalon. 

Operational  dlffleultlaa  accompanying  angina 
caaponant  off-daalgn  performance  and  throttla 
excursions  can  ba  alnlaixad  through  an  under¬ 
standing  of  tha  undarlylng  phyaleal  prlnclplaa  and 
characteristic*  of  tha  aacountarad  anomaly.  An 
understanding  of  tha  driving  aachaaiaaa  behind 
unstable  angina  oparatlon  can  proapt  hardware 
aodlf lcatlona  to  reduce  or  allalnaca  tha  source  of 
the  lnatablUty.  Coaponant  coupling  can  ralnforca 
unstable  operation  of  several  eoapoaeats  of  tha 
angina.  An  understanding  of  tha  coupling  aechanlaa 
can  result  In  design  codifications  to  alleviate 
reinforcement  of  the  Instability.  Dynaale  Insta¬ 
bilities  and  tlaa  delays  of  tha  angina  control 
eyscra  can  ba  identified  and  alnlaixad  during 
transient  operation.  An  understanding  of  tha 
phenoaena  underlying  unstable  angina  operation  la 
facilitated  through  parfotalag  appropriate  propul¬ 
sion  ayataa  flight  testa  and  applying  aethodleal 
post-flight  data  analysis. 

Propulsion  ayataa  evaluation  and  tasting  of  high 
parforaanca  aircraft  are  conducted  at  tha  Air  Pores 
Plight  Teat  Center.  During  racant  propulsion  systaa 
flight  easts,  dynaale  angina  Instabilities  have 
bean  observed  to  result  from  throttla  axcuralona 
which  Involve  transient  angaaator  oparatlon.  Those 
Instabilities  vara  recorded  by  pressure  transducers 
located  In  tha  angina  cars  and  bypass  flowpaths,  and 
In  the  gas  generator  and  augaantor  fuel  delivery 
syseaas.  Maaauraaaats  from  typical  augaantor 
transients  wars  chosen  for  analysis  cod  for  tha 
gana ration  of  a  data  base  of  tha  dynaale  angina 
Instabilities.  These  Instabilities  ware  charac¬ 
terised  by  engine  stall  or  stagnation,  augaantor 
rissblo  or  blowout,  and  pros  sure  oscillations  In 
tha  fuel  flow  aatarlng  and  distribution  systaa. 
Orla  ting  data  analysla  procedures  wars  errand  ad 
end  new  techniques  ware  developed  for  tha  analysis 
of  angina  Instabilities  racordad  la  tha  data  baaa. 
Procedures  wars  defined  which  aided  tha  Identifi¬ 
cation  and  characterisation  of  pbaooaaaa  typical 
of  unstable  augaantor  transient  oparatlon.  The 
analysis  tachnlquas  rev  salad  basic  aachaaisas 
undarlylng  tha  dynaale  Instability  as  wall  aa  a 
coupling  between  angina  conponaata.  Sir Over a 
design  aodlf lcatlona  war*  dacanalnad  Ca  allavlata 
specific  dynaale*  of  tha  fuel  delivery  ayataa  and 
improve  transient  augaantor  oparatlon. 


Propulsion  Syati  Description 

Tha  data  aaalyxad  ware  obtained  fren  a  propulsion 
systaa  which  la  typical  of  a  aodarn  high  parforaanca 
augmented  turbofan  angina.  Figure  1  Illustrates  a 
cross  sectional  achaaatle  of  this  type  of  propulsion 
systaa.  Primary  angina  stations  and  eenponancx 
are  Indicated.  Tha  angina  la  typical  of  a  low 
bypass  ratio  and  high-pressure  ratio  design.  Tha 
aultlpla-stage  fan  la  driven  by  a  aultlpla— stage, 
low-pressure  turblaa.  Similarly,  tha  compressor  la 
powered  by  a  aultlpla-scage,  high-pressure  turbine. 
Variable  compressor  vaoas  are  utilised  to  help 
aalatala  high  parforaanca  ovar  a  wide  range  of 
power  tattings.  Tha  angina  core  combustor  la 
typically  of  the  annular  type.  Tha  angina  la  also 
equipped  with  a  full  length  annular  fan  duct  and  a 


nixed  flow  augaantor.  The  augaentor  It  terminated 
by  a  variable-area,  convergent-divergent  noxsle. 
Tha  eugmencor  contains  aultlpla  spray  manifolds 
or  sa plants  which  deliver  fuel  to  tho  afterburner 
combustion  area  for  afterburner  oparatlon.  That* 
manifolds  era  located  in  tha  core  and  bypass  flow- 
paths  behind  the  low  praaaurt  turblna. 


consmko* 


STATION  numa 

Pig.  1  Croat  sectional  view  of  typical  propulsion 
systaa  with  coaponant  and  station  Identifi¬ 
cation. 

Modulation  at  tha  augaantor  fuel  delivery  syacna 
la  typically  accomplished  by  an  Integrated  after¬ 
burner  and  notxla  control  ayataa.  Tha  afterburner 
control  ayataa  has  tha  primary  function  of  governing 
tha  fuel  delivered  to  tha  afterburner  manifolds, 
and  aalatalalag  a  correct  prsasura  level  In  tha 
afterburner  by  varying  tha  noxxla  throat  eras.  Tha 
control  aytaa  la  also  responsible  for  sequentially 
supplying  the  fual  to  the  various  segaeacs  of  tha 
sugneacor.  This  prncasa  normally  coaelsts  of 
Initially  filling  tha  aaalfolds  at  a  high  rata  of 
fual  flow  (qulekflll)  and  chan  transferring  the 
aanlfnld  fuel  flow  to  a  precisely  aactrad  rata  of 
flow.  Tha  augaantor  fual  control  Is  particularly 
Important  during  throttla  excursions  which  Involve 
transient  augaantor  oparatlon. 


Data  Acquisition 

Signals  from  pressure,  taaporature,  rotation, 
and  position  sensors  wars  racordad  during  propul¬ 
sion  ayataa  flight  teats.  Tha  tnatruaencaclon  In¬ 
stalled  on  tha  aupiantad  eurbofan  angina  utilised 
for  flight  teat  la  shown  in  figure  Z.  Tha  various 
parameters  vara  aeasurad  at  tha  locations  indicated 
la  figure  Z.  Data  wars  racordad  by  a  low  frequency 
response  pulse  coda  aodulaclon  (PCM)  syacaa  sod  s 
high  frequency  response  frequency  aodulaclon  (FM) 
ayataa.  The**  data  screams  wars  racordad  oo 
magnetic  taps  aboard  tha  aireraft  during  flight 
cast  for  poat-fllght  proeaaalng.  Low  frequency 
response  paraaatara  such  s*  the  power  lever  angle 
(PLA),  augaantor  llght-off  detector  (LCD),  percent 
nossla  area  (AJ),  and  the  rotational  spaed  of  the 
high  and  low  pressure  rocors,  wars  racordad  by  the 
PCM  ayataa.  81gh  frequency  response  paraaettrs, 
such  a*  tha  core  and  bypass  duet  flovpach  pressures, 
oala  pusp  fual  pressure,  and  the  tugaaocor  sagaant 
fual  praasuraa,  ware  racordad  by  the  FM  systaa. 
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The  various  parameters  were  filtered  by  the  PCM 
a  ad  FM  lnscrusencaclon  systems  before  digitization 
to  prevenc  aliasing  trrocs. 


Fig.  2  Propulsion  syscem  Instrumentation. 


Throctls  snaps  f ran  ldls  co  maxlmun  power  ware 
performed  ac  relatively  high  aldcudes  and  low  Mach 
numbers.  Figure  4  shows  a  portion  of  the  aircraft 
flight  envelope  and  che  augmentor  operating  enve¬ 
lope.  The  augsencor  operating  envelope  la  divided 
Into  three  regions  which  are  based  on  che  probabil¬ 
ity  of  afterburner  transient  success.  Unrestricted 
afcerburner  operadon  Is  permitted  In  region  1  and 
all  augaencor  transients  are  expected  co  be  success¬ 
ful.  However,  afterburner  mlslighcs,  llght-off 
stalls  or  blowouts,  and  stalls  related  to  ruable, 
may  occur  during  throttle  transients.  Augmentor 
cranaianta  performed  In  region  3  are  likely  to  be 
unsuccessful.  Also,  afcerburner  rvseble  and  blowout 
are  possible  during  steady-state  operation.  Figure 
4  shows  a  region  of  the  augmencor  operating  envelope 
where  combustion  Instabilities  such  as  rumble  have 
been  experienced.  This  region  Is  Indicated  by  the 
cross-hatched  area. 


Data  from  che  aircraft  and  engine  lnscrioenca- 
don  recorded  on  che  alrcrafc  cape  were  processed 
by  a  series  of  dlglcal  computer  programs  which 
produced  che  fllghc  data  file.  This  process  Is 
llluscraced  In  Figure  3.  The  analog  alrcrafc 
cape  was  converted  lnco  a  digital  cape  by  an 
analog-co-dlgltal  converter.  During  this  process 
the  FM  parameters  were  digitized  ac  a  race  of 
1,000  samples  per  second.  The  data  were  then 
converted  lnco  engineering  units  and  various 
computations  were  performed  with  direction  from 
the  project  history  file.  Ac  this  time  Instru¬ 
mentation  errors  were  corrected,  and  derived  param¬ 
eters,  such  as  che  engine  pressure  rado  and 
compressor  pressure  rado,  were  calculated.  This 
processing  was  conducted  by  the  Air  Force  Flight 
Test  Cancer's  Standard  Analog  co  Digital  System 
(SANDS)  for  data  analysis  which  produced  a  digital 
data  cape  wlch  the  parameters  expressed  In 
engineering  units.  Further  processing,  such  as 
data  smoothing  and  editing,  was  performed  by  che 
Uniform  Fllghc  Test  Analysis  System  (UFTAS)  ac  che 
Fllghc  Test  Center.  UFTAS  produced  a  data  base 
compatible  wlch  dlglcal  data  analysis  software. 
The  dlglcal  daca  base  was  subsequently  referenced 
by  dlglcal  computer  programs  which  performed  analy¬ 
sis  casxa  to  Identify  and  characterize  che  measured 
augaencor  transient  phenomena. 


Fig.  3  Alrcrafc  and  post-flight  daca  acquisition 
syscem. 


Flight  tesc  results  were  obtained  from  chrottle 
snaps  from  idle  to  maximize  power  ac  operating  points 
In  regions  1  and  2  and  near  che  interface  between 
regions  2  and  3.  Idle  co  maximum  power  snaps 
were  generally  successful  In  region  1.  However, 
seme  throttle  snaps  attempted  in  region  2  led  co 
engine  stall  or  augmencor  blowout.  Throttle  snaps 
to  maximum  powar  in  region  3  frequently  resulted 
in  similar  operational  anomalies.  From  unsuccessful 
augmentor  cransiencs,  typical  flight  tesc  daca 
were  chosen  to  exemplify  transient  phenomena 
and  develop  daca  analysis  techniques.  Two  such 
flight  test  poincs  are  indicated  in  Figure  4.  A 
throctls  snap  was  attempted  at  flight  case  point  1. 
This  augmencor  transient  resulted  ia  an  engine 
stall.  A  similar  throccle  snap  was  attempted  at 
fllghc  tesc  point  2.  This  augaencor  transient 
resulted  in  an  augmencor  blowout.  Daca  obtained 
from  these  transient  events  were  recorded  and 
processed  through  the  pose-flight  data  processing 
syscem  for  further  analysis. 


Data  Analysis  Procedure 

Decs  analysis  consisted  of  processing  che  param¬ 
eter  time  hiscories  of  che  dlglcal  daca  file  through 
a  group  of  dlglcal  compueer  programs.  Figure  S 
illustrates  a  block  diagram  of  the  daca  analysis 


software.  Thia  software  la  organized  Into  four 
group*  of  functionally  alallar  progress.  Thaaa 
group*  ara  cha  data  consolidation  and  editing 
prograa  group,  the  data  dlaplay  and  Hat  pro- 
gran  group,  tho  frequency  anaiyal*  prograa  group, 
and  tho  digital  filter  prograa  group,  Th*  data 
analyala  aoftvara  utilized  tha  0PTAJ  data  file 
atruetura  for  tha  cnaannl cation  of  data  batvaan 
tha  Individual  analyala  prograaa,  Also,  fonattad 
output' files  vara  aaployad  to  tranaalt  analyala 
raaulta  to  hard  copy  davlcaa. 


rig.  J  Block  dlagran  of  data  analyala  aoftvara, 

Th*  varloua  analyala  aoftvara  groupa  contain 
prograaa  which  par  fora  apadflc  taaka  of  data 
reduction,  Th*  data  conaolldatlon  and  editing 
group  of  prograaa  porforaed  operation*  auch  aa  th* 
conaoUdntlou  of  parameter  tlaa  hlatorlea,  the  cor¬ 
rection  of  Inatruaontacloo  and  dl^ltlzln,  arrora, 
th*  ganaratlon  of  derived  paraaecara,  and  foraat 
conversion  of  data  fllea,  Parma* ter  tla*  hlatorlea 
could  be  plotted  or  Hated  la  a  aped  fled  foraat 
through  the  ueage  of  eoftvar*  la  th*  dlaplay  and 
list  prograa  group,  Th*  frequency  analyala  prograa 
group  provided  the  capability  to  analyze  data  la 
the  frequency  doaaln,  A  hat  Fourier  Transform 
(  FPT)  of  a  dyoaalc  paraaacer  or  Cha  frequency  doaaln 
tranafer  function  beevaen  tvo  paraaecara  eould  be 
dace rained.  Tha  aoftvar*  of  the  digital  flltar 
prograa  group  facilitated  tha  daalgn  and  la  plea  an¬ 
te  cl  on  of  lov-paaa  and  a  ul  tip la- band,  flnlte- 
tapulae-reeponae  digital  filter*,  A  daacrlpclon 
of  cha  digital  filter  vaa  repreeantad  by  It*  lapulae 
reponaa  which  vaa  acorad  on  th*  flltar  daacrlpclon 
flla,  Th*  aoftvar*  prograaa  within  th*  four 
prograa  group*  generated  anaiyal*  raaulta  la  the 
fora  of  llatlnga,  plot*,  or  output  fllea  which 
could  be  directed  to  a  plotter  or  11a*  printer. 


Oeta  Analyala  Result* 


technique*  that  are  dlacuaaed  Include  POt  and  PM 
data  else  hletory  plot*,  Fourier  Tranatom  based 
apectral  analyala,  and  digital  filtering.  Thaa* 
analyala  tool*  aaalated  la  revealing  th*  charsc- 
tarlacle*  and  phyalcal  baa la  of  dyoaalc  Instabil¬ 
ities  of  aupientor  operation  that  were  observed 
during  flight  teat. 


PCM  Data  Tla*  Hlatorlea 

Tha  tie*  histories  of  angina  paraaecara  recorded 
by  th*  PCM  dace  acquisition  eyacea  provided  a 
caaprahenalve  record  of  a  transient  event  during 
fHght  cast.  Aa  advantage  of  thia  data  acquisition 
syataa  1*  that  aany  paraaecara  eould  be  Input  to  tha 
POt  data  screaa.  Therefore,  an  Interaction  or 
coupUng  beevaen  engine  components  could  be  dis¬ 
cerned  during  analysis,  PCM  daea  of  eh*  engine 
paraaaeara  were  recorded  during  eh*  augneneor  tran¬ 
sient  flight  testa,  figures  6  and  7  present  PCM 
daea  obtained  froa  flight  ease  points  1  and  2. 
Various  paraaecara  are  presented  over  a  10  second 
tlaa  Interval  which  capeures  the  transient  processes 
loading  to  a  stall  or  blowout. 
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Pig.  4  Tlaa  history  of  PCM  data  obtained  during 
au^eator  transient  at  flighe  case  point  1. 

Tha  transient  processes  which  occurred  during  th* 
throttle  snap  to  aaxlatas  power  at  flight  teet 
point  l  resulted  In  an  engine  stall.  These  tran¬ 
sient  processes  are  characterized  by  th*  paraneter 
tlaa  hlatorlea  given  la  Figure  6.  Th*  PLA  was 
advanced  to  the  naxlaia  power  sectlng  at  the 
elapsed  tla*  of  t  ■  0  seconds.  Approximately  8 
seconds  later,  sfeer  full  thrust  augmentation  <** 
achieved,  so  engine  stall  was  experienced  and 
tha  PLA  was  cutback  to  an  lncemedlata  power 
sectlng. 


Results  froa  various  data  analysis  technique*  Tha  paraaacer  tla*  histories  presented  In  Figure 
chat  were  eppUed  to  augnontor  transient  data  6  sueoarlz*  cha  sequence  of  evencs  which  occurred 
are  presented  and  dlacuaaed  below.  Oata  analysis  between  the  chrotcl*  snap  and  the  Incipient  stall. 


Varlou*  praaaurt*  la  cha  fan  and  cora  duel,  ai  wall 
a  a  cha  axhauac  nozzle  parcaac  area  variation,  ara 
shown  along  with  cha  augaancor  fuel  manifold  prat- 
suras*  Ac  c  «  0  saconds  cha  PLA  requasc  raaulcad 
la  an  lneraasa  of  cha  fan  duet  prasaura  (Pt2.5c), 
cora  flow  prasaura  bahiad  cha  high  prassurs  coa- 
praaaor  (Pb),  and  prasaura  la  cha  mixed  flow  raglon 
bahiad  cha  low  prasaura  Curblna  (Pt6m).  Tha  rlsa 
la  cha  aagina  flowpaeh  prasauraa  gradually  irr- 
craasad  uadi  c  •  3.3  saconds.  Ac  chit  ciae  cha 
augaancor  f  ual  coocrol  bagan  co  sequentially 
praaaurlza  cha  augaancor  fual  spray  rings  which 
cooaiscad  of  sagaaac  1  through  sagaanc  3.  Augaen- 
cor  lgnldon  ac  t  •  4  saconds  was  daducad  fron  cha 
rapid  lneraasa  In  Pc6a  and  confirmed  by  cha  LCD. 
Afcar  sagaanc  1  was  scabla,  cha  fual  control  qulck- 
flllad  sagaanc  2  ac  c  •  3.3  saconds  and  chan  trans- 
farrad  sagaanc  2  co  aacarad  flow.  This  procass  was 
concinuad  for  sagsanes  3  chrough  3  and  camlnacad 
with  sagaanc  3  on  aacarad  flow  ac  approximately 
t  •  7.7  saconds.  Tha  cora  and  fan  flowpaeh  pres¬ 
sures  Increased  and  decreased  In  aapllcuda  during 
chls  procass  In  response  co  backpressure  genaraced 
by  Cha  augaancor  scappad  coabuadon  procassas  and 
modulation  of  cha  axhauac  nozzle.  Pc6a  and  Pt2.5c 
experienced  substantial  variations  In  aapllcuda 
during  cha  sequential  accivaeion  of  cha'  augaancor 
segments.  However,  any  Increase  In  backpressure 
Co  cha  fan  up  to  C  *  7  saconds  was  accommodated  by 
cha  fan' s  stall  margin.  Also,  cha  nozzle  area 
adjusted  co  maintain  cha  corracc  engine  prasaura 
rado  ( EPR)  •  Afcar  cha  quickflll  and  transfer  co 
macared  flow  of  sagaanc  3,  a  prasaura  pulse  was 
generated  which  had  sufficient  aapllcuda  co  lneraasa 
tha  fan  pressure  ratio  co  a  value  above  lea  sceady- 
staca  scan  line.  This  raaulcad  In  a  fan  stall  ac 
t  •  8  saconds-  Tha  stall  was  evidenced  by  a 
rapid  lneraasa  In  Pt6m  and  a  savara  oscillation 
of  Pc2.3c.  immediately  following  cha  stall,  cha 
PLA  was  broughc  Co  an  Intermediate  power  setclng, 
cha  fual  flow  was  sequentially  stopped  co  cha  aug- 
mancor  segaancs,  cha  aozzia  was  closed,  and  cha 
angina  recovarad  ac  c  ■  9.3  saconds. 

An  augaancor  blowouc  was  also  axparlancsd  during 
a  chrocda  excursion  from  Idle  co  maximal  power. 
The  blowouc  occurred  ac  flight  case  point  2.  Figure 
7  summarizes  the  time  historian  of  cha  PCM  param¬ 
eters  during  10  saconds  of  time  which  capcured  cha 
evenc.  Tha  flva  augaancor  sagaanc  fual  prassuras 
are  shown  along  wlch  cha  preseures  measured  In  cha 
fao  and  cora  flowpaeh  and  cha  parcaac  exhaust 
aozzia  area  variation. 

Ac  c  •  0.8  seconds  of  elapaad  time,  cha  chrocda 
was  snapped  from  an  idle  power  setting  co  cha  max¬ 
imise  power  setting  of  cha  PLA.  Again,  cha  chrocda 
advance  resulted  In  an  lneraasa  In  cha  aapllcuda 
of  cha  cora  and  fan  flowpaeh  prassuras.  Ac  c  *  4 
seconds  cha  augaancor  sagiants  wars  sequentially 
prasaurlzad.  Ignition  of  Cha  augaancor  was  detect¬ 
ed  ac  c  >4.3  saconds.  Sagaancs  2  chrough  3 
proceeded  chrough  a  sequence  of  qulckfills  and 
craasfers  co  macared  flow  uncll  e  •  8.3  seconds 
whan  sagaanc  3  was  supplied  wlch  aacarad  fual  flow. 

Subscandal  variations  In  cha  aapllcuda  of  cha 
fan  duct  prassuras  ware  Induced  due  co  prasaura 
pulses  from  cha  augaancor  cransienc  fual  flow. 
However,  cha  stall  margin  of  cha  fan  was  not 
axcaadad  and  cha  nozzle  adjuscad  co  maintain  a 
constant  angina  prasaura  ratio.  One  second  afcar 


sagaanc  3  vs a  transferred  to  macarad  flow,  ac 
c  -  9.3  seconds,  an  augaancor  blowouc  was  exper¬ 
ienced.  The  causa  of  chls  blowouc  was  difficult 
co  determine  from  Cha  PCM  data  alone* 


i >U  AUGNENTON 

AOVANCE  IGNITION  SLOWOUT 


0123«SS  7  3  3*  10 


EUMED  TINE  (SECONDS) 


Fig.  7  Tima  history  of  PCM  daca  obcainad  during 
augaancor  transient  ac  fllghc  test  point  2. 

Tha  PCM  daca  acquisition  system  provided  a 
comprehensive  summary  of  the  variation  of  many  an¬ 
gina  paramacars  during  augaancor  transients.  Tha 
Interaction  between  angina  components  was  described 
by  a  large  number  of  paramacars  over  a  long  period 
of  clma.  However,  the  lowar  sample  races  associated 
with  cha  PCM  daca  syscem  llmlcad  cha  ability  Co 
Identify  rapid  fluctuations  In  paramecer  values. 
From  cha  PCM  daca,  lc  was  difficult  co  dacataine  cha 
phenomena  which  Induced  tha  pressure  oscillacions 
Chat  raaulcad  In  an  angina  scall  and  blowouc.  A 
better  undarscandlng  of  the  augaancor  transient 
evancs  was  obcainad  chrough  the  FM  time  history 
daca. 


FM  Data  Tima  Histories 

Tha  FM  daca  acquisition  syscaa  provided  paramecer 
time  histories  which  could  ba  digitized  ac  a  high 
race.  A  digitizing  race  of  1,300  samples  per  second 
was  employed.  Characteristics  of  several  augaancor 
cransienc  phenomena  in  cha  gas  flowpaeh  and  fual 
delivery  syscem  ware  revealed  ac  Chls  high  sample 
race. 

The  FM  Oaca  revealed  decails  of  augaancor  crao- 
slanc  phenomena  which  pracadad  the  englna  scall 
char,  occurred  ac  flight  cesc  point  1.  The  time  his¬ 
tories  of  several  FM  parameters  ara  given  in  Figure 
8  for  cha  scall  event.  These  time  histories  span 
a  period  of  1  second  in  time  which  corresponds  to  an 
elapaad  time  of  tha  PCM  daca  presented  in  Figure  6 
from  c  -  7.2  seconds  to  c  *  3.2  saconds.  Tha 
quickflll  of  tha  augaancor  segaancs,  as  exemplified 


by  (•(■•at  5  «c  c  •  7.4  seconds,  1(  shown  co  con¬ 
sist  of  high  amplitude  pressure  fluctuations 
which  dlalnlshad  la  saplltud*  whan  th«  sagaant 

trsasfarrad  to  a  aatarad  flow  rata.  Tha  dynoalcs 
ssaoclatad  with  cha  qulekflll  of  tha  au^ancor 
aaalfolds  war*  transalttad  to  sagaancs  already 
supplied  with  aatarad  fuel  flow  dua  co  cha  coanoa 
delivery  sysesa.  This  Is  raalizad  by  eoaparlag 
cha  qulekflll  of  ••(sane  5  and  cha  rasulcsoc 
pressure  osclllstloos  In  tha  aatarad  fuel  flow 
of  sagsaot  3.  Also,  cha  transfer  to  aatarad 

fuel  flow  was  not  saooch  sad  was  rrraalad  to 

ba  characterized  by  a  prtssurs  splka.  This  prasaura 
splka  has  bean  reported  to  craata  a  praaaura 

pulsa  In  tha  fan  duct  flowpath  of  sufflelant  am- 
plituda  co  Indue*  a  fan  stall.  Tha  sagaanta  already 
on  aacarad  flow  also  azparlaacad  a  praaaura  drop 
slaca  they  are  supplied  by  cha  coonon  fuel  daliv- 
•  ry  systaai. 
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Fig.  8  Tiaa  history  of  FM  data  pracadlng  angina 
stall  at  flight  Cast  point  1. 

Tha  qulekflll  and  tranafar  to  aacarad  flow  of 
cha  augaancor  sagaants  war*  obaarvad  to  Indue* 
prasaura  oaclllaclona  In  cha  angina  fan  duct  aa 
shown  by  Flgurs  8.  This  phanoaanon  Is  axampliflad 
by  tha  varladon  of  cha  duct  staclc  prasaura  at 
station  5,  P*5c,  during  tha  qulekflll  of  sagaant  3. 
Tha  induced  prasaura  oaclllaclona  In  tha  fan  flow- 
path  bacaaa  aora  savara  as  tha  qulekflll  process 
was  applied  by  subaaquant  aanlfold  sagaants. 
This  la  dua  co  tha  common  fuel  dalivary  systaa 
of  cha  aaalfolds.  Tha  induced  gas  flowpath  prasaurt 
oaclllacioaa  war*  obaarvad  with  reduced  saplltud* 
by  other  flowpath  pressure*  aaasurad  la  cha  fan 
duct  and  also  by  Pt2.3h  which  was  aaasurad  la  tha 
core  flowpath  behind  cha  fan.  Sowavar,  no  Influence 
was  obaarvad  In  n  aaaaurad  downstream  of  cha 
hlgh-prassura  compressor.  This  suggaats  chat  tha 
flowpath  prasaura  oscillation*  war*  noe  crsnsmlc- 
ted  through  tha  conprassor. 

Figure  8  also  shows  tha  onset  of  praaaura 
oscillations  In  tha  fan  duct  at  t  •  7.7  seconds 
which  tiaaadiaealy  pracadad  stall.  7hasa  prasaura 


oscillations  raachad  an  amplitude  which  caused 
cha  fan  co  stall  at  t  ■  8  seconds.  During  this 
period  of  claa,  tha  augaancor  fuel  dalivary  systaa 
appeared  to  ba  stabilized  with  all  sagaants  on 
aacarad  flow.  Thasa  pressure  oscillations  In  cha 
gas  flowpath  ware  obaarvad  through  cha  fan  duct 
and  also  at  tha  hlgh-prassura  compressor  Inlet- 
Howavar,  cha  duet  flow  pressure  oscillations  war* 
noe  obaarvad  at  tha  axle  of  cha  high-pressure 
coaprassor  as  shown  by  Pb.  Again,  cha  pressure 
oscillations  war*  obaarvad  to  ba  readily  transalt¬ 
tad  through  tha  fan  but  not  through  tha  high-pres- 
■ura  compressor.  This  Information  suggaats 
that  tha  large  saplltud*  prasaura  fluctuations 
aay  ba  associated  with  a  combustion  instability 
of  tha  augaantor  which  occurred  after  sagaant  3 
stabilized  with  aatarad  fuel  flow.  Combustion 
Instabilities  of  this  type  can  ba  drlvan  by  a 
coupling  baevaan  an  acoustic  oscillation  at  a 
resonant  node  of  tha  augaancor  and  tha  hast  addi¬ 
tion  process. 

Tha  transient  phenomena  preceding  cha  augsencor 
blowout  axparlancad  at  flight  test  point  2  are 
eharaetarlzad  by  tha  FM  parameters  presenetd  In 
Flgura  9.  Thasa  data  correspond  to  evo  periods  of 
time  recorded  by  Cha  PCM  data  system  and  plotted  In 
figure  7.  Thasa  claa  periods  correspond  co  (lapsed 
class  from  c  •  8.2  seconds  to  c  -  8.7  seconds  and 
from  t  ■  9.2  seconds  to  t  -  9.7  seconds.  Tha  data 
prasancad  In  Figure  9  from  c  -  8.2  seconds  co  t  1 
3.7  second*  daaonstrata  tha  efface  of  a  modifica¬ 
tion  to  cha  augaancor  qulekflll  systaa  which 
reducad  tha  amplitude  of  prasaura  oscillations  In 
Cha  fual  dalivary  system.  Tha  time  history  of 
Psag3  shows  a  raducad  prasaura  oscillation  dua  to 
qulekflll  and  a  smoother  transfer  to  aacsred  flow. 
Also,  tha  variation  of  Psag3  indicates  chat  pres¬ 
sure  oscillations  its  no  longer  readily  trans¬ 
ferred  through  cha  fual  dalivary  system  of  tha 
augaancor. 
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Fig.  9  Tima  history  of  FM  data  preceding  augmancor 
blowout  ac  flight  test  point  2. 


The  faa  duct  pressure  oscillations  which  pre¬ 
ceded  augaeator  blowout  are  also  shown  in  Figure  9 
during  the  period  of  elapsed  time  frors  t  *  9.2 
seconds  to  t  ■  9.7  secoads-  During  this  period  of 
time,  the  augaentor  fuel  delivery  systes  had  stabil¬ 
ized  and  all  cha  segments  were  on  metered  flow.  How¬ 
ever,  pressure  oscillations  were  observed  la  the 
fan  duct  and  ac  the  coapressor  face  with  diminished 
amplitude.  As  with  flight  test  point  1,  these  tiae 
history  plots  suggest  that  a  coabustion  instability 
of  che  augaentor  had  induced  Che  blowout  during 
che  throttle  transient. 

the  above  augaentor  transient  tiae  histories 
have  suggested  several  phenomena  which  Involve 
oscillacioas  of  che  augaentor  fuel  delivery  syscea 
and  Induced  oscillations  of  che  fan  and  core 
flowpach.  These  phenomena  appear  very  similar, 
occur  within  short  time  Intervals,  and  result  in 
cwo  different  operational  anomalies :  an  engine 
stall  and  an  augaeator  blowout.  It  is  not  certain 
by  viewing  che  data  time  histories  if  the  pressure 
oscillations  of  the  fuel  system  and  flowpach  result 
from  the  same  phenomenon,  cwo  coupled  phenomena, 
or  separate  mechanisms.  A  quantitative  method 
for  analyzing  che  time  history  data  was  needed. 
Che  approach  was  to  analyze  the  digitized  EM  data 
in  che  frequency  domain. 


Power  Spectral  Density  Analysis 

Some  characteristics  of  the  mechanisms  under¬ 
lying  dynamic  data  are  noc  easily  determined  in 
Che  time  domain.  It  can  ba  more  appropriate  to 
transform  che  parameter  time  history  loco  the 
frequency  domain  and  analyze  the  signal.  The 
digitized  EM  paraaecer  time  histories  of  flight 
test  point  1  were  transformed  into  the  frequency 
domain  and  power  spectral  densities  of  che  data 
were  obealned.  These  power  spectral  densities 
identified  component  frequencies  of  the  signal 
oscillations  and  che  relative  amplitude  of  che 
peculations  at  various  engine  locations. 

Power  spectral  density  plots  of  PsSc  were  ob¬ 
tained  during  che  augaentor  quickflll  and  imme¬ 
diately  preceding  stall.  Figure  10  shows  chat  che 
characteristic  frequency  of  the  fuel  pressure  oscll- 
lacion  during  quickflll  differs  from  cne  character¬ 
istic  frequency  of  the  fan  flowpach  dynamics  prior 
to  scall.  The  fuel  pressure  dynamics  are  shown 
to  be  characterized  by  pressure  oscillations  at 
approximately  45  Hertz.  This  phenomenon  was 
observed  throughout  the  augmencor  fuel  delivery  sys¬ 
tem,  transmitted  through  the  augaentor  combustion 
process  and  resulted  in  pressure  oscillations  in 

the  fan  duct.  The  fuel  sysem*  dynamics  associated 
with  che  quickflll  process  were  reduced  in  amplitude 
after  segment  5  transferred  to  metered  fuel  flow 
as  shown  in  Figure  8.  Immediately  following 
segment  5  transfer  co  metered  flow.  Figure  10 
indicates  chat  pressure  oscillacioas,  sc  a 
characteristic  frequency  of  60  Hertz  and  lesser 
amplitude  compared  co  quickflll,  timsediacely 
preceded  che  fan  scall.  It  is  suspeccad  chat 
an  augsancor  combustion  lnscablllcy,  which  is 
characterized  by  longitudinal  mode  frequencies 
from  50  co  100  Hertz,  Induced  this  fan  stall. 
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Fig.  10  Power  spectral  density  of  Ps5c  during 
augaentor  quickflll  and  pre-stall  dynamics 
at  flight  test  point  l. 

Augaentor  ruable  was  further  suspected  since  the 
amplitudes  of  the  pressure  oeclllatloas  were  great¬ 
est  ac  station  5  in  che  fan  duct.  Figure  11  illue- 
scraces  a  comparison  of  power  spectral  density 
plots  obtained  from  measurements  ac  various  engine 
scaclons.  The  pressure  osclllselon  measured  by  Ps5c 
had  reduced  in  amplitude  ac  the  fan  exit  plane  as 
shown  by  Pt2.Sc.  The  amplitude  of  the  pressure 
oscillation  was  reduced  further  when  measured  in 
che  core  flowpach  at  scadon  2.5.  Also,  the  ampli¬ 
tude  of  the  spectral  density  plot  of  ?b  shows  chat 
the  pressure  oscillation  was  noc  readily  transmitted 
through  the  compressor. 
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Fig.  11  Power  spectral  density  of  fan  and  core 
flowpach  parameters  preceding  engine  stall 
at  flight  cast  poinc  1. 

A  pach  of  propagation  is  suggested  by  the  ampli¬ 
tude  of  the  power  speccral  densities  at  che 
characteristic  frequency  of  che  combuscion  insta¬ 
bility.  The  source  of  the  preseure  oeclllation  due 
to  tvnble  was  in  che  augmencor.  This  pressure  oscil¬ 
lation  propagated  upscream  In  che  fan  duct.  Upon 
reaching  che  fan  face,  the  pressure  oscillation 
propagaetd  further  upstream  and  was  also  traasmic- 
Cfd  through  the  faa  core  flowpach.  However,  che 
pressure  oscillation  was  absorbed  by  che  high- 
pressure  compressor.  This  resulted  in  stable 
pressure  levels  ac  che  core  burner  inlet. 


Digital  rilt«r  Processing 

Tha  fraquancy  spectrvaa  of  a  signal  often  repre¬ 
sents  cha  characteristics  of  tha  signal  aora  expll- 
cltly  than  lta  claa  domain  counterpart.  Tha  ran  oval 
or  enhancement  of  signal  eharaetarlatlea  la  tha 
frequency  domain  eaa  halp  Identify  tha  phanoaana 
under  Investigation.  An  example  of  this  la  tha 
reaoval  of  nolsa  due  to  a  power  supply  from  other- 
vise  elaan  and  daelrabla  data.  By  filtering  tha 
band  of  freqaaelaa  that  con-ala  tha  nolaa  or  a 
s taller  undesirable  band  of  frequencies*  a  aora 
useful  signal  eao  be  raaeorad.  This  signal  ean  than 
ba  analysed  without  confusion  or  lntarfaranea  (roa 
Irrelevant  frequency  eontaat. 


to  aspect  5  of  cha  augaantor.  Tha  pre-stall  dy- 
oaalca  observed  In  flowpach  praaaurea  are  shown  to 
Figure  13.  During  this  period  of  elan*  PsegS  was 
stable  and  was  being  supplied  by  aacerad  flow. 
These  pre-stall  dyuaaies  ara  probably  assoelaeed 
with  an  augaantor  eoabuatloa  instability  such  as 
nabla.  Flowpach  pressure  oaelllaclons  ware 
aaasured  at  station  i  In  the  tan  duct  (PsSc)  with 
tha  largest  aaplltuda.  As  the  pressure  oaellla- 
tlona  propagated  up  cha  fan  duct  to  the  fan  face 
and  downacraaa  through  the  fan,  they  were  aaasured 
by  Pt2.5h  la  the  angina  core  flowpach  at  a  reduced 
aaplltuda.  These  oscillations  were  not  substan¬ 
tially  tranaalcted  through  cha  high  pressure  com¬ 
pressor*  as  ft  Indicates.  This  analysis  enhanced 
cha  signals  of  Interest. 


Digital  decs  filtering  tachnlquas  wars  applied 
co  cha  paraaecer  claa  histories  obtained  during  cha 
augaantor  craaaleoc  operation.  Application  of  tha 
digital  filter  co  tha  claa  history  data  required 
soaa  knowledge  of  tha  characteristics  of  tha  data  as 
wall  as  cha  specifications  of  cha  filter  to  ba 
eaployad.  The  dynamic  paraaaters  aaasured  during 
tha  aupiaacor  transients,  ware  analysed  In  cha 
range  of  frequencies  of  Interest.  This  range  of 
frequencies  was  determined  by  cha  power  spectral 
density  of  tha  parameter.  Tha  digital  filter 
passband  was  specified  co  paaa  a  band  of  frequen¬ 
cies  containing  tha  characteristic  fraquancy  of 
the  fuel  syscaa  qulckflll  praaaura  oaelllaclons, 
pre-stall  flowpach  praaaura  dynamics,  and  eha 
dynamics  of  cha  flowpach  presaurae  before  blowout. 
A  3-pole  luccervorth  filter  uaa  designed  and 
laplMantad  with  a  forward  and  reverse  paaa  through 
cha  data  co  cancel  any  phase  shift  due  to  tha 
filter. 

*<■•  dynamic  parameters  required  conditioning 
With  a  masking  technique  before  being  fllcered. 
This  reduced  any  Induced  ringing  of  cha  filter 
from  rapid  data  drop-off.  Taramacar  HJe  obtained 
sc  flight  test  point  1  exhibited  this  rnpld  drop-off 
of  amplitude  Immediately  following  Cha  dynamics  of 
interest.  Figure  12  illustrates  cha  application 
of  a  masking  signal  co  parameter  Pa 5c  before  being 
filtered.  A  aasklng  signal  (MA SC)  was  ganaraced 
and  subtracted  from  tha  original  signal  (ORIG) 
to  produce  an  oscillating  type  of  signal  (AC) 
which  was  free  of  data  drop-off  for  procaaslng. 
Tha  aasklng  signal  consists  of  a  constant  ampli¬ 
tude  signal  uncll  tha  point  of  data  drop-off  whare 
a  copy  of  tha  original  signal  la  used  co  termin- 
ace  cha  signal.  The  resultant  AC  signal  was 
filtered  into  cha  various  bands  of  interest  without 
error  due  to  filter  ringing. 

The  digital  filtering  technique  described  above 
was  applied  co  cha  dynamic  parameters  measured  at 
flight  cast  point  1.  Figure  13  Illustrates  cha 
digital  filter  procaaslng  results  for  tha  flowpach 
paraaaters  PsSe,  Ft2.3h,  Fb,  and  tha  augaantor  fuel 
manifold  pressure  of  segment  3.  These  paraaatars 
ware  filtered  In  a  fraquancy  band  from  14  Harts  co 
70  Berts  and  over  a  period  of  time  which  snearn  passed 
cha  qulckflll  dynamics  of  Cha  fuel  system  and  the 
pre-stall  flowpath  dynamics.  This  period  of  clma 
corresponds  to  an  elapsed  time  la  Figure  i  from  c  • 
7.2  to  c  •  8.2.  Tha  pressura  oaelllaclons  in  PsegS 
dua  co  Cha  qulckflll  process  ara  clearly  shown  and 
are  observed  to  Influence  cha  fan  duct  pressure, 
FsSe,  end  the  core  flowpach  pressure,  Pt2.5h.  This 
influence  was  probably  dua  to  a  caabuaclon  Insta¬ 
bility  driven  by  oscillating  fuel  delivery  rates 
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Fig.  12  ficaapl*  of  digital  fllttr  ttchtilqua. 
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Fig.  13  Dlglcal  filter  results  of  flowpach  and  fuel 
syscaa  parameters  at  flighc  cesc  point  l. 


A  similar  digital  (Hearing  analysis  was  conducc- 
ed  on  prsssura  parameters  measured  In  cha  augsantor 
fual  da 11 vary  syscem  and  gas  flowpath  prior  to  cha 
blowout  ancouncarad  ac  (light  ease  point  2.  Xha 
raaules  shown  In  Figure  14  correspond  to  an  elapsed 
time  (roa  c  -  8.5  seconds  to  c  *  9.5  seconds. 
Ourlng  this  period,  the  augsantor  (ual  prasaura  was 
stable.  However,  as  Flgura  14  illustrates  substan¬ 
tial  pressure  oscillations  wars  axparlancad  In  cha 
gas  (lowpaeh.  Paraaatars  Ps5c  and  Pt2.5h  axparl¬ 
ancad  ehasa  oscillations  and  chair  ralaelve  magnl- 
cudaa  suggest  ehac  Che  source  was  an  augaancor 
combustion  Instability.  This  relationship  is 
easily  deceralned  through  cha  enhanced  signals 
produced  by  cha  digital  (Hearing  procedure. 
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Fig.  L4  Olgical  (liter  results  o(  (lowpaeh  and  (uel 
syscaa  paraaatars  obtained  at  (light  ease 
point  2. 

Olgical  (Hearing  techniques  ware  ea cloyed  to 
determine  if  cha  characteristic  frequency  o (  the 
augaancor  Instability  was  constant  (or  a  period  oi 
ciaa  iomadlacaly  preceding  stall.  Ps5c  was 
analyzed  during  an  elapsed  ciaa  of  c  -  7.2  seconds 
to  c  •  8.2  seconds  for  cha  augaancor  cranslanc 
axparlancad  ac  flight  case  poinc  1.  Flgura  15 
shows  cha  results  from  filtering  PsSc  into  four 
frequency  bands  which  range  from  30  Hertz  to  70 
Hertz.  Ourlng  cha  initial  growth  of  cha  pressure 
oscillation,  tha  higher  frequency  bands  con¬ 
tained  moat  of  cha  frequency  concent.  However, 
as  tha  stall  point  approached,  Che  characteristic 
frequency  of  cha  (lowpaeh  oscillation  was  shifted 
to  the  lower  frequency  bands  and  was  aost  pro¬ 
nounced  jusc  before  stall.  Tha  shift  in  eharac- 
carlsclc  frequency  was  also  revealed  by  a 
frequency  analysis  of  cha  ducc  pressure  fluct¬ 
uations  lseaadlacsly  preceding  stall-  Flgura  16 
Illustrates  this  frequency  shift  pracadlng  stall. 
Tha  characteristic  frequancy  is  seen  co  shift  to 
a  lower  frequency  just  prior  to  scall. 
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Fig.  15  Multiple  frequency  band  digital  filter 
results  of  Ps5c  during  pre-stall  dynamics 
at  flight  test  point  1. 


Fig.  16  Variation  in  characteristic  frequency  of 
pre-stall  dynamics  ac  flight  test  point  1. 

Tha  shift  in  characteristic  frequency  occurred 
when  Che  augaancor  was  suspected  co  be  experiencing 
r<sbla.  This  frequency  shift  say  indicate  chat  cha 
frequency  of  the  combustion  Instability  was  ap¬ 
proaching  a  resonant  sods  of  the  augaancor  system. 
Jusc  prior  co  scall,  this  resonant  soda  coupled 
with  an  acoustic  osclllacion  co  produce  large 
pressure  oeclllacions,  which  subeaquencly  stalled 
cha  fan. 


Concluding  Remarks 

A  method  has  been  developed  for  the  analysis  of 
gas  turbine  engine  dynamic  instabilities  chac  have 
been  obaerved  during  transient  augaancor  operation. 
These  instabilities  were  characterized  by  engine 
scall  or  stagnation,  augaancor  ratable  or  blowout, 
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aad  praaaura  oscillations  to  cho  augaantor  fuel 
distribution  ayacaa.  Gristing  data  analysis 
procedures  vara  axtaodad  aod  new  taehalquaa  vara 
developed  for  eha  analysis  of  these  angina 
Instabilities  which  vara  racordad  by  praaaura 
aanaora  1 oca cad  to  tha  aogloa  flowpaeh  and  fual 
delivery  ayacaa.  Analysis  procedures,  locludad 
paraaacar  tlaa  hlacory  aaalyala  of  law  aod  high 
frequency  raaponaa  data,  Fourier  Traoaf ora  baaad 
apaccral  aaalyala,  aod  data  flltarlag.  Thaaa 
aaalyala  taehalquaa  revealed  aschanlsas  underlying 
the  dynaalc  lnatabllielaa  aad  loearacelooa  batwaeo 
angina  components. 

a  caaprahanalva  hlatory  of  augoantor  traoaleot 
avanea  vaa  suaaarlzed  by  low  fraquaocy  raaponaa 
data  obtained  through  eha  PCI  data  aequlaleloo 
systaa.  Many  paraaetara  ware  Input  to  the  PCM  data 
aeraaa  which  captured  tha  overall  eraoalaoe  process— 
aa  and  coaponanc  loearacelooa  pracadlog  a  fao  atall 
or  aupMotor  blowout.  However,  dua  eo  eha  low 
data  saapla  rata,  deeaila  of  tha  obaervad  aooualy 
ware  not  obealoed  by  eha  PCM  data  aequlaleloo 
syatao. 

The  PM  data  aequlaleloo  eyataa  provided  paraoatar 
else  hiatorlaa  which  could  ba  digitized  at  a  high 
aaopla  rata.  Oecallad  paraoatar  tlaa  hlatory  analy- 
ala  haa  auggaatad  aavaral  phaoaoaoa  which  la- 
volva  oaclllatlona  of  eha  augoantor  fual  delivery 
syatao  and  the  augaaotor  eoabuaeloo  procaea.  Tha 
augment  aanlfold  qulchflll  procaaa  waa  determined 
eo  Induce  fual  praaaura  oaclllatlooa  which  ware 
eranaateted  through  eha  fual  delivery  ayaean  eo 
ocher  aanlfold  assents.  Also,  thaaa  fual  praaaura 
oacillaelooa  induced  flowpaeh  praaaura  oscillations 
through  the  augoantor  caabuacioa  procaaa.  Au^an- 
tor  enable  waa  indicated  through  largo  amplitude 
praaaura  oaclllatlona  In  tha  angina  core  and  fao 
flowpaeh. 

Analysis  of  cha  digitized  IH  parameter  tlaa 
hiatorlaa  In  eha  frequency  doaoln  provided  a  aachod 
of  dlacarnlng  tha  charactarlaclc  frequency  of  eha 
augaantor  operating  anaoaly.  ulffarant  characear- 
latle  fraquandaa,  for  the  fual  ayacaa  dynnalca  and 
eoabuaeloo  lnatabllity  pracoding  stall  or  blowout, 
suggaaead  that  tha  fual  ayacaa  dynaalca  . ware  not 
directly  related  to  cha  augoantor  caabuacioa 
dynaalca.  a  propagation  path  of  tha  core  aad  fan 
duct  praaaura  oaclllatlona,  which  originated  In 
tha  augaantor,  waa  deteralaed  by  the  amplitude  of 
the  paraaacar  power  apaccral  denalclee  at  eha 
characearlacic  fraquaocy  of  eha  eoabuaeloo  lnaea- 
blllcy. 

iaaoval  or  anhancaraant  of  elgnal  charzctarlaclca 
In  cha  frequency  doaala  waa  parforaad  by  utilizing 
digital  filtering  cechnlquaa.  Tha  enhanced  param¬ 
eter  tlaa  hiatorlaa  revealed  tha  Influence  of  fual 
syscaa  praaaura  oaclllatiooa  on  eha  gaa  praaauras 
In  tha  fan  aad  cora  flowpaeha.  Augoantor  combus¬ 
tion  praaaura  oaclllatlooa  ware  shown  eo  be  easily 
cranaolctad  through  eha  fan  but  not  through  tha 
compressor.  Tha  characteristic  frequency  of  the 
augaantor  combustion  Instability  waa  shown  eo  ba 
dm a  dependant.  Tha  fraquaocy  shift  prior  to 
stall  Indicated  a  coupling  of  a  resonant  aoda  of 
tha  coabuatloa  lnatabllity  with  an  acoustle 
oscillation  eo  produce  large  pressure  oaclllatlona 
to  tha  fan  duct. 


CONVERSION  OF  A  B-FILE  TO  A  C-FILE 


PURPOSE: 

BTOC  CONVERTS  AN  UFTAS  DATA  FILE  FROM  A  B-FILE  FORMAT  TO  A  C-FILE 
FORMAT.  THE  PROGRAM  ENABLES  THE  SPECIFICATION  OF  THE  C-FILE  HEADER 
RECORD  INFORMATION. 

REQUIRED  INPUT: 

BTOC  CAN  BE  RUN  BOTH  INTERACTIVELY  AND  BATCH.  INTERACTIVE  EXECUTION 
IS  DOCUMENTED  HERE.  A  LOCAL  B-FILE  NAMED  BFILE  IS  THE  PRIMARY  INPUT 
FOR  BTOC.  THIS  FILE  MUST  CONTAIN  THE  PARAMETER  TIME  HISTORIES  IN  A 
B-FILE  FORMAT  FOR  THE  CONVERSION  TO  A  C-FILE  FORMAT.  CONTROL  OF  PRO¬ 
GRAM  EXECUTION  IS  OBTAINED  THROUGH  AN  EXECUTION  CARD  STATEMENT  WHICH 
IS  ENTERED  AT  THE  TERMINAL  AS  A  COMMAND.  THE  FORMAT  OF  THIS  COMMAND 
IS: 

BTOC 

OPTIONAL  INPUT: 

A  LOCAL  FILE  NAMED  INPUT  CAN  BE  USED  TO  SPECIFY  INFORMATION  TO  BE 
INCLUDED  ON  THE  OUTPUT  C-FILE  HEADER  RECORD.  THIS  INFORMATION  IS 
READ  BY  BTOC  THROUGH  THE  FOLLOWING  PROGRAM  RECORDS: 

READ  (5,5,END=10)  LI C , IT AI L , I TEST , FLT , DFLT , DR EQ , DCOM 
5  FORMAT  (A5,2X,I3.2XtI1,2X,A5,2X,A6,2X,A6,2X,A6) 

READ  (5. ' (10(A6) ) ' ,END=20)  (TITLE(I) ,1*1 , 10) 

OUTPUT: 

BTOC  PRODUCES  A  C-FILE  NAMED  CFILE  WHICH  CONTAINS  THE  INPUT  PARA¬ 
METER  TIME  HISTORIES  IN  A  C-FILE  FORMAT. 


CONSOLIDATION  OF  C-FILE  PARAMETERS 


PURPOSE: 

CONSO  CONSOLIDATES  TIME  HISTORY  PARAMETERS  FROM  MULTIPLE  C-FILE 
B-FILE  SECTIONS  AND  OUTPUTS  THESE  PARAMETERS  TO  ONE  B-FILE  OR 
C-FILE.  THIS  CONSOLIDATION  OCCURS  OVER  A  SELECTED  TIME  PERIOD  WHICH 
MUST  BE  COMMON  TO  ALL  B-FILE  SECTIONS  OF  THE  INPUT  C-FILES.  MODERATE 
DATA  DROPOUT  DUE  TO  DATA  COMPRESSION  IS  AUTOMATICALLY  CORRECTED. 

REQUIRED  INPUT: 

1.  NUMBER  OF  C-FILES  TO  BE  CONSOLIDATED. 

2.  START  TIME-OF-DAY  FOR  THE  RESULTANT  B-FILE  (SECONDS). 

3.  FINISH  TIME-OF-DAY  FOR  THE  RESULTANT  B-FILE  (SECONDS). 

4.  DATA  SAMPLE  RATE  (SAMPLES/SECOND). 

5.  NAMES  OF  THE  INPUT  C-FILES  TO  BE  CONSOLIDATED. 

6.  SECTION  NUMBER  OF  THE  B-FILES  TO  BE  CONSOLIDATED. 

7.  NAME  FOR  THE  OUTPUT  B-FILE  OR  C-FILE. 

OUTPUT: 

1.  B-FILE  OR  C-FILE  WHICH  CONTAINS  THE  CONSOLIDATED  PARAMETERS. 


MANIPULATION  AND  CONDITIONING  OF  PARAMETER  SEQUENCES 


PURPOSE: 

SEQCOND  MANIPULATES  PARAMETER  TIME  HISTORIES  THAT  ARE  INPUT  FROM  A 
C-FILE  AND  RECOVERS  MODERATE  DATA  DROPOUT  DUE  TO  COMPRESSION  FROM 
SANDS  PROCESSING.  DATA  SEQUENCE  CONDITIONING  CAPABILITIES  INCLUDE: 
AMPLIFICATION  OR  ATTENUATION  OF  THE  SEQUENCE  AND  ITS  DC  LEVEL,  THE 
DERIVATION  OF  NEW  PARAMETERS,  AND  THE  APPLICATION  OF  A  COSINE  TAPER 
TO  THE  END  OF  A  SEQUENCE.  THE  RESULTANT  PARAMETER  TIME  HISTORIES  ARE 
OUTPUT  TO  A  SPECIFIED  LOCAL  B-FILE. 

REQUIRED  INPUT: 

1.  NAME  OF  THE  INPUT  C-FILE. 

2.  B-FILE  SECTION  NUMBER  CONTAINING  THE  PARAMETERS  TO  BE  MANIPULA¬ 
TED. 

3.  NUMBER  OF  PARAMETERS  TO  BE  CONDITIONED. 

4.  NUMBER  OF  PARAMETERS  TO  BE  DERIVED. 

5.  DATA  SAMPLE  RATE  (SAMPLES/SECOND). 

6.  START  TIME-OF-DAY  FOR  THE  RESULTANT  TIME  HISTORIES  (SECONDS). 

7.  FINISH  TIME-OF-DAY  FOR  THE  RESULTANT  TIME  HISTORIES  (SECONDS). 

8.  NAME  FOR  THE  B-FILE  WHICH  RECEIVES  THE  CONDITIONED  DATA. 

9.  FOR  EACH  THROUGHPUT  PARAMETER  TO  BE  CONDITIONED: 

A.  PARAMETER  NAME. 

B.  DC  BIAS  (IF  DESIRED). 

C.  ATTENUATION  FACTOR  (IF  DESIRED). 

D.  COSINE  TAPER  (IF  DESIRED). 

10.  FOR  EACH  DERIVED  PARAMETER: 

A.  PARAMETER  NAME. 

B.  CHOICE  OF  THE  DERIVED  PARAMETER  TO  BE  A  RATIO  OR  DIFFERENCE 
OF  PARAMETERS. 

OPTIONAL  INPUT: 

1.  SELECTION  OF  A  COSINE  TAPER  REQUIRES: 

A.  TIME-OF-DAY  OF  THE  SIGNAL  APEX  (SECONDS). 

8.  APPROXIMATE  PERIOD  OF  THE  SIGNAL  DYNAMICS  (SECONDS). 


OUTPUT: 


1.  LOCAL  B-FILE  WHICH  CONTAINS  THE  MANIPULATED  DATA  SEQUENCES  AS  A 
TRANSFERRED  OR  DERIVED  PARAMETER. 


TIME  SHIFT  OF  A  B-FILE  PARAMETER  SEQUENCE 


PURPOSE: 

TIMALIN  PRODUCES  A  B-FILE  TIME  HISTORY  WHICH  IS  SHIFTED  IN  TIME  REL¬ 
ATIVE  TO  THE  ORIGINAL  TIME  HISTORY.  THE  SHIFTED  PARAMETER  SEQUENCE 
CAN  BE  APPENDED  TO  THE  INPUT  B-FILE  AS  AN  ADDITIONAL  PARAMETER  OR 
WROTE  OVER  THE  INPUT  PARAMETER.  THE  INPUT  SIGNAL  MUST  BE  BAND  LIMIT¬ 
ED  AND  NOT  WINDOWED. 

REQUIRED  INPUT: 

1.  NAME  OF  THE  LOCAL  B-FILE  WHICH  CONTAINS  THE  PARAMETER  TO  BE  TIME 
SHIFTED. 

2.  NAME  OF  THE  PARAMETER  TO  BE  SHIFTED. 

3.  START  TIME-OF-DAY  FOR  THE  PARAMETER  TIME  HISTORY  (SECONDS). 

4.  FINISH  TIME-OF-DAY  FOR  THE  PARAMETER  TIME  HISTORY  (SECONDS). 

5.  DATA  SAMPLE  RATE  (SAMPLES/SECOND). 

6.  DESIRED  TIME  SHIFT  (SECONDS). 

7.  NAME  FOR  THE  LOCAL  B-FILE  WHICH  CONTAINS  THE  RESULTANT  PARAMETER 
TIME  HISTORY. 

8.  SELECTION  OF  THE  RESULTANT  TIME  HISTORY  TO  BE  A  NEW  OR  OLD  PARA¬ 
METER. 

OPTIONAL  INPUT: 

1.  SELECTION  OF  THE  RESULTANT  TIME  HISTORY  TO  BE  A  NEW  PARAMETER 
REQUIRES: 

A.  NAME  FOR  THIS  NEW  PARAMETER. 

OUTPUT: 

1.  B-FILE  WHICH  CONTAINS  THE  TIME  ALIGNED  PARAMETER  TIME  HISTORY  AS 
A  NEW  OR  OVERWRITTEN  PARAMETER. 
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B-FILE  OR  C-FILE  LISTING  OF  PARAMETERS 


PURPOSE: 


LISTF  IS  A  MODIFIED  VERSION  OF  AN  UFTAS  UTILITY  PROGRAM  WHICH  PRINTS 
THE  CONTENTS  OF  A  B-FILE  OR  C-FILE.  THE  PROGRAM  ENABLES  A  SELECTION 
OF  THE  DATA  PRINT  RATE,  PARAMETERS  TO  BE  LISTED,  AND  THE  DATA  FOR¬ 
MAT. 


REQUIRED  INPUT: 

LISTF  CAN  BE  RUN  BOTH  INTERACTIVELY  AND  BATCH.  INTERACTIVE  EXECUTION 
IS  DOCUMENTED  HERE.  CONTROL  OF  PROGRAM  EXECUTION  IS  OBTAINED  THROUGH 
AN  EXECUTION  CARD  WHICH  IS  ENTERED  AT  THE  TERMINAL  AS  A  COMMAND.  THE 
FORMAT  OF  THIS  COMMAND  IS: 

LISTF (FILE) 

WHERE  FILE  IS  THE  NAME  OF  A  LOCAL  FILE  WHICH  CONTAINS  THE  PARAMETERS 
TO  BE  LISTED.  THE  DEFAULT  NAME  IS  CFILE. 

OPTIONAL  INPUT: 

LISTING  OPTIONS  ARE  INVOKED  THROUGH  AN  ARGUMENT  LIST  ON  THE  EXECU¬ 
TION  CARD.  THIS  IS  SPECIFIED  AS: 


LISTF(FILE,PR=N1,PC=N2,SD=N3) 

WHERE: 

1.  PR  IS  THE  PRINT  RATE.  N1  IS  A  NUMERICAL  VALUE  SPECIFYING  EVERY 
N1TH  POINT  TO  BE  LISTED.  DEFAULT  IS  PR=1. 

2.  PC  IS  A  FLAG  FOR  PARAMETER  SPECIFICATION.  IF  N2=0,  ALL  PARAMETERS 
ARE  PRINTED.  IF  N2=1,  A  LIST  OF  PARAMETERS  IN  FREE  FORMAT  IS  READ 
FROM  A  LOCAL  FILE  NAMED  INPUT.  ONLY  THESE  PARAMETERS  ARE  LISTED. 
DEFAULT  IS  PC=0. 

3.  SD  IS  A  FLAG  FOR  THE  NUMBER  OF  DATA  VALUE  SIGNIFICANT  DIGITS.  IF 
SD=0,  5  SIGNIFICANT  DIGITS  ARE  LISTED.  IF  SD=1,  7  SIGNIFICANT 
DIGITS  ARE  LISTED.  DEFAULT  IS  SD=0. 


OUTPUT: 


THE  LISTING  OF  DATA  VALUES  IS  WRITTEN  TO  A  LOCAL  FILE  NAMED  OUTPUT. 


PLOT  OF  A  B-FILE  PARAMETER  TIME  HISTORY 


PURPOSE: 

PLOT  DISPLAYS  A  PLOT  OF  A  SPECIFIED  PARAMETER  TIME  HISTORY  ON  THE 
SCREEN  OF  A  TEKTRONIX  GRAPHIC  TERMINAL. 

REQUIRED  INPUT: 

1.  MODEL  NUMBER  OF  THE  TEKTRONIX  GRAPHIC  TERMINAL. 

2.  NAME  OF  THE  B-FILE  WHICH  CONTAINS  THE  DATA  SEQUENCE  TO  BE 
PLOTTED. 

3.  NAME  OF  THE  B-FILE  PARAMETER  TO  BE  PLOTTED. 

4.  UNITS  OF  THE  B-FILE  PARAMETER. 

5.  DATA  SAMPLE  RATE  OF  THE  SEQUENCE  (SAMPLES/SECOND). 

6.  PLOT  START  TIME-OF-DAY  (SECONDS). 

7.  PLOT  FINISH  TIME-OF-DAY  (SECONDS). 

OPTIONAL  INPUT: 

1.  STATEMENT  OF  A  PLOT  TITLE  OR  REMARK. 

2.  SELECTION  OF  MANUAL  PLOT  SCALING  REQUIRES: 

A.  PLOT  ORDINATE  MAXIMUN  VALUE. 

B.  PLOT  ORDINATE  MINIMUM  VALUE. 


OUTPUT: 

1.  PLOT  OF  THE  B-FILE  PARAMETER  TIME  HISTORY  ON  THE  SCREEN  OF  A 
TEKTRONIX  GRAPHIC  TERMINAL. 


FAST  FOURIER  TRANSFORM  OF  A  B-FILE  PARAMETER 


PURPOSE: 

FFT  PERFORMS  A  FAST  FOURIER  TRANSFORM  ON  AN  INPUT  B-FILE  PARAMETER 
TIME  HISTORY.  THE  FFT  RESULTS  CAN  BE  LISTED  AT  THE  TERMINAL  SCREEN 
WRITTEN  TO  AN  OUTPUT  FILE,  OR  PLOTTED  ON  THE  SCREEN  OF  A  TEKTRONIX 
GRAPHIC  TERMINAL. 

REQUIRED  INPUT: 

1.  MODEL  NUMBER  OF  THE  TEKTRONIX  TERMINAL  (IF  APPROPRIATE). 

2.  NAME  OF  THE  B-FILE  WHICH  CONTAINS  THE  DATA  SEQUENCE  TO  BE  AN¬ 
ALYZED. 

3.  NAME  OF  THE  PARAMETER  TO  BE  ANALYZED. 

4.  UNITS  OF  THE  PARAMETER. 

5.  DATA  SAMPLE  RATE  (SAMPLES/SECOND). 

6.  ANALYSIS  START  TIME-OF-DAY  (SECONDS). 

7.  NUMBER  OF  POINTS  TO  BE  ANALYZED. 

8.  SELECTION  OF  THE  FORMAT  FOR  THE  FFT  RESULTS. 

OPTIONAL  INPUT: 

1.  APPLICATION  OF  A  HANN  WINDOW  TO  THE  DATA  PRIOR  TO  ANALYSIS. 

2.  SELECTION  OF  OUTPUT  FORMAT  REQUIRES: 

A.  CHOICE  OF  TERMINAL  SCREEN  OR  OUTPUT  DUMP  FILE. 

3.  SELECTION  OF  GRAPHIC  DISPLAY  OF  RESULTS  REQUIRES: 

A.  SELECTION  OF  THE  TYPE  OF  PLOT  (MAGNITUDE, PHASE, OR  PSD  VS. 
FREQUENCY). 

B.  SELECTION  OF  AUTOMATIC  OR  MANUAL  SCALING. 

4.  STATEMENT  OF  PLOT  A  REMARK. 

OUTPUT: 

1.  TABULATION  OF  FFT  RESULTS  ON  THE  TERMINAL  SCREEN  OR  ON  A  SPECI¬ 
FIED  OUTPUT  FILE. 

2.  PLOT  OF  FFT  RESULTS  AT  THE  TERMINAL  SCREEN. 


TRANSFER  FUNCTION  BETWEEN  TWO  B-FILE  PARAMETERS 


PURPOSE: 

TRANS  COMPUTES  THE  TRANSFER  FUNCTION  BETWEEN  TWO  B-FILE  PARAMETER 
TIME  HISTORIES.  RESULTS  CAN  BE  LISTED  AT  THE  TERMINAL  SCREEN,  DUMPED 
TO  A  LOCAL  FILE,  OR  PLOTTED  AT  THE  SCREEN  OF  A  TEKTRONIX  GRAPHIC 
TERMINAL. 

REQUIRED  INPUT: 

1.  MODEL  NUMBER  OF  THE  TEKTRONIX  GRAPHIC  TERMINAL  (IF  APPROPRIATE). 

2.  NAME  OF  THE  LOCAL  B-FILE  WHICH  CONTAINS  THE  PARAMETERS  TO  BE  ANA¬ 
LYZED. 

3.  ANALYSIS  START  TIME  OF  DAY  (SECONDS). 

4.  NUMBER  OF  POINTS  FOR  THE  ANALYSIS. 

5.  DATA  SAMPLE  RATE  (SAMPLES/SECOND). 

6.  NAMES  OF  THE  PARAMETERS  TO  BE  ANALYZED. 

7.  SELECTION  OF  AN  OUTPUT  FORMAT  FOR  THE  ANALYSIS  RESULTS. 

OPTIONAL  INPUT: 

1.  CHOICE  OF  RESULTS  LISTED  AT  THE  TERMINAL  SCREEN  OR  DUMPED  TO  A 
LOCAL  OUTPUT  FILE. 

2.  SELECTION  OF  GRAPHIC  DISPLAY  OF  RESULTS  REQUIRES: 

A.  TYPE  OF  PLOT  (MAGNITUDE  OR  PHASE  VS.  FREQUENCY). 

B.  CHOICE  OF  AUTOMATIC  OR  MANUAL  SCALING  OF  THE  PLOTS. 


OUTPUT: 

1.  TABULATION  OF  TRANS  RESULTS  AT  THE  TERMINAL  SCREEN  OR  ON  A  SPEC¬ 
IFIED  OUTPUT  FILE. 

2.  PLOTS  OF  TRANS  RESULTS  ON  THE  SCREEN  OF  A  TEKTRONIX  GRAPHIC  TER- 


DESIGN  OF  LOWPASS  DIGITAL  FILTERS 


PURPOSE: 

LOWPASS  IS  A  DESIGN  TOOL  FOR  THE  GENERATION  OF  MAXIMALLY  FLAT  PASS- 
BAND  AND  STOPBAND  FINITE  IMPULSE  RESPONSE  (FIR)  DIGITAL  LOWPASS  FIL¬ 
TERS.  THE  FILTER  DESIGN  IS  USER  SPECIFIED  AND  THE  RESULTANT  FILTER 
DESCRIPTION  IS  WRITTEN  TO  A  FILTER  FILE  IN  STANDARD  FORMAT.  THIS 
FILTER  DESIGN  HAS  A  LINEAR  PHASE  VARIATION  WITH  FREQUENCY  WHICH  IM¬ 
PARTS  A  CONSTANT  TIME  DELAY  TO  THE  FILTERED  DATA. 

REQUIRED  INPUT: 

1.  SAMPLE  RATE  OF  THE  DATA  TO  BE  FILTERED  (SAMPLES/SECOND). 

2.  CENTER  FREQUENCY  OF  THE  TRANSITION  BAND  (HZ). 

3.  WIDTH  OF  THE  TRANSITION  BAND  (HZ). 

4.  NAME  FOR  THE  FILTER  DESCRIPTION  FILE. 

5.  SELECTION  OF  LISTING  OPTIONS  FOR  THE  FILTER  COEFFICIENTS. 

OPTIONAL  INPUT: 

1.  SELECTION  OF  A  LISTING  OF  THE  FILTER  COEFFICIENTS  REQUIRES: 

A.  CHOICE  OF  A  TERMINAL  LISTING  OR  DUMP  TO  A  SPECIFIED  LOCAL 
FILE. 

OUTPUT: 

1.  LOCAL  FILTER  FILE  FOR  SUBSEQUENT  APPLICATION  TO  DATA  THROUGH 
THE  FILTER  IMPLEMENTATION  PROGRAM. 

2.  IF  CHOSEN,  A  LOCAL  FILE  WHICH  CONTAINS  A  LISTING  OF  THE  FILTER 
COEFFICIENTS  AND  THE  FILTER  DESCRIPTION. 


DESIGN  OF  MULTIPLE  PASSBAND  A  STOPBAND  DIGITAL  FILTERS 


PURPOSE: 

MLTPASS  IS  A  DESIGN  TOOL  FOR  THE  GENERATION  OF  EQUIRIPPLE  MULTIPLE 
PASSBAND  AND  STOPBAND  FINITE  IMPULSE  RESPONSE  (FIR)  DIGITAL  FILTERS. 
THE  FILTER  DESIGN  IS  USER  SPECIFIED  AND  THE  RESULTING  FILTER  DES¬ 
CRIPTION  IS  WRITTEN  TO  A  SPECIFIED  LOCAL  FILE  IN  STANDARD  FORMAT. 

THE  RESULTANT  FILTER  DESIGN  HAS  A  LINEAR  PHASE  VARIATION  WITH  FRE¬ 
QUENCY.  THIS  IMPARTS  A  CONSTANT  TIME  DELAY  TO  THE  FILTERED  DATA. 

REQUIRED  INPUT: 

1.  FILTER  LENGTH. 

2.  TOTAL  NUMBER  OF  PASSBANDS  AND  STOPBANDS. 

3.  FILTER  GRID  DENSITY. 

4.  SAMPLE  RATE  OF  THE  DATA  TO  BE  FILTERED  (SAMPLES/SECOND). 

5.  SPECIFICATIONS  FOR  EACH  PASS3AND  AND  STOPBAND: 

A.  BAND  LOWER  EDGE  FREQUENCY  (HZ). 

B.  3 AND  UPPER  EDGE  FREQUENCY  (HZ). 

C.  GAIN  OF  THE  BAND. 

D.  RELATIVE  WEIGHT  OF  THE  BAND. 

6.  NAME  FOR  THE  FILTER  DESCRIPTION  FILE. 

7.  SELECTION  FROM  LISTING  OPTIONS  FOR  A  DESCRIPTION  OF  THE  FILTER 
SPECIFICATIONS. 

OPTIONAL  INPUT: 

1.  SELECTION  OF  THE  FILTER  DESCRIPTION  TO  3E  WRITTEN  TO  A  FILTER 
FILE  REQUIRES: 

A.  NAME  FOR  THE  FILTER  FILE. 

2.  SELECTION  OF  A  FILTER  DESCRIPTION  LISTING  REQUIRES: 

A.  CHOICE  OF  A  TERMINAL  LISTING  OR  LISTING  ON  A  LOCAL  FILE. 

OUTPUT: 

1.  LOCAL  FILTER  FILE  WHICH  CONTAINS  THE  FILTER  DESCRIPTION  AS  RE¬ 
QUIRED  BY  THE  FILTER  IMPLEMENTATION  PROGRAM. 

2.  IF  CHOSEN,  A  LOCAL  OUTPUT  FILE  WHICH  CONTAINS  A  LISTING  OF  THE 
FILTER  DESCRIPTION  AND  THE  FILTER  COEFFICIENTS. 


DIGITAL  FILTER  IMPLEMENTATION 


PURPOSE: 

FILTER  PERFORMS  THE  FILTERING  OF  DIGITAL  DATA  BY  EMPLOYING  A  FINITE 
IMPULSE  RESPONSE  (FIR)  FILTER  DESIGN  WHICH  HAS  BEEN  PREVIOUSLY  GEN¬ 
ERATED.  A  SPECIFIED  B-FILE  PARAMETER  IS  FILTERED  AND  APPENDED  TO 
THE  INPUT  B-FILE.  THIS  FILE  IS  OUTPUT  AS  A  UNIQUE  B-FILE. 

REQUIRED  INPUT: 

1.  NAME  OF  THE  B-FILE  WHICH  CONTAINS  THE  PARAMETER  TO  BE  FILTERED. 

2.  NAME  FOR  THE  OUTPUT  B-FILE. 

3.  NAME  OF  THE  PARAMETER  TO  BE  FILTERED. 

4.  NAME  FOR  THE  FILTERED  PARAMETER. 

5.  SELECTION  OF  THE  FILTER  DESCRIPTION  TO  COME  FROM  A  FILTER  FILE 
OR  TERMINAL  ENTRY. 

OPTIONAL  INPUT: 

1.  SELECTION  OF  THE  FILTER  DESCRIPTION  FROM  A  FILTER  FILE  REQUIRES: 
A.  NAME  OF  THE  FILTER  FILE. 

2.  SELECTION  OF  THE  FILTER  DESCRIPTION  FROM  THE  TERMINAL  REQUIRES: 

A.  LENGTH  OF  THE  FILTER. 

B.  FILTER  GAIN  FOR  EACH  SEQUENCE  POINT. 

OUTPUT: 

1.  B-FILE  CONTAINING  THE  FILTERED  PARAMETER  AS  AN  APPENDED  PARAMETER 
TO  THE  INPUT  B-FILE. 


WHITE  NOISE  DATA  FILE 


PURPOSE: 

WHITE  IS  A  WHITE  NOISE  DATA  B-FILE  WHICH  PROVIDES  A  SOURCE  TEST  SIG¬ 
NAL  WHICH  IS  CHARACTERIZED  BY  FREQUENCY  CONTENT  UNIFORMLY  SPREAD 
ACROSS  THE  SAMPLED  DATA  SPECTRUM.  THIS  TIME  HISTORY  CAN  BE  USED  TO 
TEST  A  DIGITAL  FILTER  DESIGN  UNDER  CONSIDERATION.  THE  IMPLEMENTATION 
OF  THE  FILTER  TO  THIS  WHITE  NOISE  FILE  YIELDS  OUTPUT  WHICH  INDICATES 
THE  FREQUENCY  RESPONSE  OF  THE  FILTER. 

SPECIFICATIONS: 

1.  CONTAINS  ONE  PARAMETER  NAMED  NOISE1. 

2.  DATA  SAMPLE  RATE  IS  1000  SAMPLES /SECOND. 

3.  START  TIME-OF-DAY  OF  THE  SEQUENCE  IS  0.001  SECONDS. 

4.  FINISH  TIME-OF-DAY  OF  THE  SEQUENCE  IS  1.024  SECONDS. 

5.  LENGTH  OF  THE  SEQUENCE  IS  1024  SAMPLES. 

6.  NOISE1  CONTAINS  501  COMPONENT  SINUSOIDS  (INCLUDES  DC). 

7.  FREQUENCY  SPACING  BETWEEN  SINUSOIDS  IS  1.0  HZ  BETWEEN  0  AND 
500  HZ. 


